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PREFACE 


This unit shows you ways in which the computer can extend your knowledge 
and understanding of introductory quantum mechanics. The unit presents 
numerical solutions to Schrodinger’s equation; the solutions are easy to under- 
stand and yet include solutions to many Schrodinger equation problems 
which can’t be handled by analytic methods. The emphasis in the unit will 
be on ways the computer can extend the quantum mechanics you under- 
stand. The unit assumes some knowledge of introductory quantum mechanics, 
specifically analytic solutions for the infinite square well, and matching wave 
functions at boundaries. 


Your instructor can assign chapters and exercises out of this unit in much the 
same way as out of a textbook. Lectures on the material may be presented, 
but you should be able to understand the material without lectures. The 
harder exercises can be used as the basis of projects if your instructor wishes. 


Each chapter in the unit starts with a brief discussion of the physics dis- 
cussion of the physics discussed in the chapter and then moves on to an 
explanation of the numerical procedure used with the computer. Exercises 
follow with one exercise completely worked out in the text so that you can 
see what a sample solution looks like. This sample exercise is sometimes a 
problem for which the answer is already known and therefore provides an 
extra check on the computer method. Some exercises are marked with 
asterisks. One asterisk means the exercise is fairly demanding while two aster- 


isks indicates a very challenging problem. Numerical methods will not be 


discussed much in the text; there is a section in the appendix which dis- 
cusses the half-step integration used in the first two chapters. 


There are, of course, other problems concerning quantum mechanics which 
can be done with a computer. On the one hand, the computer can calculate 
energies and wave functions for analytic solutions like the finite square well 
or the harmonic oscillator or the hydrogen atom. On the other hand; the 
computer can sometimes do the matrix operations which occur in perturba- 
tion theory or in the matrix formulation of quantum mechanics. These other 
applications are not included because this unit emphasizes extensions of 
introductory quantum mechanics byond ordinary analytic methods. 


There are also ways to solve the Schrodinger equation directly as a partial 
differential equation using the computer. These applications will not be in- 
cluded because of their extra sophistication. When dealing with partial ditfer- 
ential equations on a computer, you have to worry about convergence, 
stability, and representation. Even though the one dimensional, time- 
dependent Schrodinger equation is a relatively simple partial differential 
equation of the parabolic type, the most straightforward approaches using a 
Computer turn out to be unstable. The solutions for such problems are 
discussed in reference (1). 


The three-dimensional, time-independent Schrodinger equation is an example 
of an elliptic partial differential equation (as are Laplace’s equation and 
Poisson’s equation). These solutions are more straightforward; examples of 
elliptic partial differential equations are discussed in references (2) and (3). 
Most of the interesting potentials of quantum mechanics have symetries 
which do not allow you to use partial differential equation techniques. 


Throughout the booklet the atomic units introduced by Hartree (Proc. Camb. 
Phil. Soc. 24, 89 (1928)) will be used. The units have the advantage that the 
numbers being calculated stay near 1; numbers rarely get too small or too 
large for the computer. Atomic units are often used in quantum mechanics; 
the units have the effect of setting h = e = m, = 1. In atomic units. 
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guantity unit 

charge e 

mass m, 

distance first Bohr radius (.562917A) 

velocity first Bohr speed (2.2 X 10® cm/sec) 

energy twice the ionization potential of hydrogen 
(27.21eV) 

time first Bohr period/27(2.419 X 10°'” sec) 


Often graphical output is useful. If you have a plotter or a CRT terminal 
available, fine. If not, terminal plotting (using the teletype to plot a graph) 
can be helpful. The appendix gives a program to convert printing programs 
{programs which type out lists of numbers) to programs which plot on the 
terminal. This terminal plotting will be used several times in the unit to give 
you several examples. 


John Merrill 
Hanover, New Hampshire 


This unit was written while the author was on the faculty of the Department 
of Physics and Astronomy at Dartmouth College, Hanover, New Hampshire. 
John Merrill has used computers to support his classwork and research since 
joining the Dartmouth faculty in 1966. During the period of 1970-72 he was 
part of a team of professors who developed computer classroom use in engin- 
eering, mathematics and physics under the auspices of NSF funded project 
COEXIST. 


John has written many journal articles on specific uses of computers tn 
physics teaching, particularly in the two main physics teaching journals, The 
Physics Teacher and The American Journal of Physics. He has written several 
booklets on computer use in physics teaching including another unit in the 
Hewlett-Packard series, Electric and Magnetic Fields. John is also the author 
of the book, Computers in Physics, soon to be published by Houghton- 
Mifflin Company. He is now the Director of the Center for Educational Design 
at Florida State University. 


Special credits go to Gregory Hughes, a recent PhD in physics from Dart- 
mouth who helped develop the exercises and to Chris Doerr of Hewilett- 
Packard who learned quite a bit of quantum mechanics while copy editing 
the manuscript. 





TABLE OF CONTENTS 


WrRPANMN- = 


wk 


19 


19 
20 
22 
31 


35 


CHAPTER ONE: ONE-DIMENSIONAL BOUND STATES 
Introduction 

The Computer Approach 

Initialization 

The Harmonic Oscillator 

Exercises 


CHAPTER TWO: THREE-DIMENSIONAL BOUND STATES 
Introduction 

The Equations in @ and ¢ 

The Radial Equation 

Initialization 

Exercises 


CHAPTER THREE: ONE-DIMENSIONAL, TIME-DEPENDENT 
SCHRODINGER EQUATION FOR TIME-INDEPENDENT 
POTENTIALS 

Introduction 

States of the Infinite Square Well 

The Behavior of Sets of Discrete States 

One-Dimensional Continuum States 

Free Wave Packets 

The Propagation of Free Wave Packets 

The Propagation of a Packet over a Potential Well 

Exercises 


CHAPTER FOUR: THREE-DIMENSIONAL CONTINUUM STATES 
FOR SPHERICALLY SYMMETRIC POTENTIALS 

Introduction 

Boundary Conditions for P(r) 

The Iterative Procedure 

Continuum States for a Screened Coulomb Potential 

Scattering Theory 

Exercises 


APPENDIX 

Introduction 

Terminal Plotting 

Symmetry Proof for One-Dimensional Potentials 
The Half-Step Iterative Integration 


SELECTED ANSWERS 


REFERENCES 


TABLE OF CONTENTS UC 


Vv 








ONE-DIMENSIONAL BOUND STATES UU 1 


CHAPTER ONE: ONE-DIMENSIONAL BOUND STATES 


~ INTRODUCTION 


The one-dimensional, time-independent, Schrodinger equation for the wave 
function, P, of a particle in a potential, V(x), is 


Ah? d?P(x) 
+ P E 

ea eee V(x)P(x) = EP{x) (1) 

where E is the energy of the state with wave function P. In Hartree atomic 

units (see the Preface) this equation becomes 


1 d?P 
- > 2 + V(x)P = EP (2) 


or 
p’’ = 2(V(x) -&) P (3) 


When V(x) > E for both large, positive x and large, negative x, a discrete set 
of allowed energies occurs. This set occurs because, when E > V(x), the wave 
function is oscillatory; this is the ‘‘classically allowed’’ region of space. When 
E < V(x), the wave function decays toward zero; this is the ‘‘classically for- 
bidden” region of space since, classically, when the total energy, E, of the 
particle is less than the potential energy, V(x), the kinetic energy of the 
particle must be negative. In classical physics, a negative kinetic energy, 
Yemv2, implies an imaginary velocity, v, which is the mathematical way of 
saying such motion cannot occur (in classical physics). 


In quantum mechanics, the particle can spend some time in regions of space 
where E < V(x), These are the classically forbidden regions outside the edges 
of a finite square well, for example. Quantum mechanical tunneling through 
barriers is also a result of this quantum mechanical penetration of classically 
forbidden regions of space. 


Consider a region of space where E < V(x) and where V(x) * a constant V,, 
The Schrdédinger equation (3) becomes p"’ = 2(Vo- -E)P. The two solutions 
to this equation (two because the equation is a second order differential 
equation) are Pae “OX where @ = V 2(V,-E). If the region where E < V(x) is 
bounded, both solutions can occur, but “if the region is unbounded, then one 
of the two solutions is not allowed. If the region is unbounded towards +, 
then e*®* grows without bound; this implies that IP|2 (the probability 
density for the particle) is arbitrarily large the farther out you go. Such an 
infinite probability of finding the particle at infinity is not allowed physically. 
Hence, the wave function must decay towards zero (in this case exponentially) 
as x becomes large. A similar argument holds for x tending toward ~°o; now it 
is the e“* term which diverges. The result is still that, as |x | gets large, P 
tends to zero. 


These demands on the wave function’s behavior at large positive and negative 
x are the boundary conditions used both by the analytical methods of solu- 
tion (those in which the solution is an equation for P) and by the numerical, 
computer-based methods presented in this booklet. Because the wave function 
must go smoothly to zero at large x, the particle can have only certain, dis- 
crete (i.e. separated from each other) total energies. As you will see, for other 
energies, the wave functions will not go to zero at large x. 
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P(x) 





= an eigenvalue 


P(x) 





One of the big differences between analytical methods (leaving out approxi- 
mate methods you learn in higher level quantum mechanics) and the method 
presented here is that analytical approaches can only solve five quantum me- 
chanical problems (the infinite square well, the finite square well, the finite 
square barrier, the harmonic oscillator, and the hydrogen atom). The com- 
puter method can solve almost any Schrodinger equation problem for almost 
any potential, V. In this booklet we will limit ourselves to symmetric poten- 
tials since the vast majority of interesting potentials are symmetric (i.e. are 
such that V(—x) = V(x)). 


THE COMPUTER APPROACH 


The computer (or numerical) method is an “‘algorithmic” (or “iterative” or 
“step-by-step’’) method of solution. The computer approach integrates the 
Schrodinger equation by taking small steps, Ax, and finding the new values 
of the wave function, P, and its derivative, P’ = dP/dx, after each little step. 


Suppose you know the wave function, P, and its first derivative, P’ at some 
point x. Suppose you also know how P should behave at some other x (for 
example, that P should go to zero as x gets large). Given any potential, 
V(x), you try an energy, E; you set the initial values of P and P’ at the point 
they are known; you take a sma!! step Ax and calculate the new values of P 
and P’ by using P(x+Ax) = P(x) + P’(xtAx/2)Ax and P’(x+Ax) = P’(x) + 
P’’ (x+Ax/2)Ax. Schrodinger’s equation gives you the value of the second deriva- 
tive, P’’ = d2P/dx2. You now use the same procedure again to go from x+tAx 
to the point x+2Ax, and so forth. In this way, you walk out x (integrating 
the Schrodinger equation step-by-step as you go) until you reach the other 
place at which you know the behavior of P. You then just look at the be- 
havior of P; if E was one of the allowed energies (eigenvalues), then P will 





behave correctly (e.g., P will go smoothly to zero for large x); if E was not a 
correct energy, then P will not behave (e.g., P will diverge for large x). If E 
is just more than an eigenvalue, P will diverge towards one infinity; if E is 
just below that eigenvalue, P will diverge to the other infinity. You can home 
in on eigenvalue energies quite quickly. 


Figure 1 shows a block diagram (a simplified computer flow chart) for this 
procedure. One possible BASIC implementation of this procedure is: 


1g DEF nvixo=....., potential, V(x) 

29 PRINT “END OF INTERVAL, DELTA-X» PARITY (@=EVENS 1=ODD). ENERGY?" 
3@ INPUT X9-X7sP.E 

a@ LET Pee-Pi=xo=8 

$@ LET X6=x9/26 rere 

6@ IF Pel THEN 92 Initialization 

79 LET Péel 

6G GOTO 1e¢@ 

98 LET Pi#1 


166 «=6LLET reneeceiereentmert Integration of the 


-Schrodinger equation 


118 LET P@ePGe ( Pi+P2eX7/2)4XT 

128 LET Pl#Pi+PeexX7 

138 LET X@#XO+X7 

14@ IF x@<x6 THEN 178 

158 LET X6"X64X9/26 Print loop 
168 PRINT XO. POs POePa 


176 IF x6<x9 THEN 18@— Return for next Ax step 


188 PRINT 
199 @oTo 28— Return for new parameters 
2668 END 
. Ne } ‘ 
aa : ‘ 
{ or r 
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Figure 1. 






Initialize 
PP’, x 






Calculate 
? 
new P,P, x 


Flow Chart of Computer Approach One- 
Dimensional Quantum Mechanics 
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A few comments on the program are in order: First, there is no one “‘right”’ 
implementation of the procedure. Any program that calculates the quantities 
correctly is “right.” The program shown is just one such correct implementa- 
tion. Second, notice the use of the half-step method (see the Appendix) in 
the lines calculating P2 = d2P/dx? from the Schrodinger equation and P@ = the 
wave function, P, from the first derivative, P1 = dP/dx. Third, the print-test 
lines are not necessary, but it is convenient not to print out on the termina! 
the results of all several hundred steps, Ax, you will be taking. Fourth, the 
potential, V(x), is given as a defined function, FNV(X), to emphasize that 
this short program solves any (symmetric potential) Schrodinger equation 
problem. All you do is type the potential into the DEF line. Finally, the 
initialization of P@ = the wave function and P1 = the first derivative of the 
wave function with respect to x will be discussed in the next section. 


INITIALIZATION 


Much of upper level quantum mechanics hinges on the development of 
theorems concerning the Schrodinger equation (or its equivalents). One such 
theorem states: 


The wave functions for the stationary states (i.e. the solutions of the 
time-independent Schrodinger equation) for a symmetric potential, 


V(x), can always be written as completely even (i.e., P(—x} = P(x)) or 
completely odd (i.e., P(-—x) = —P{x)). 





At first glance that theorem (which is proved in the Appendix for those of 
you with sufficient mathematical background) has little to do with our prob- 
lem. In fact, the theorem gives our initial values of P and P’ atx =0. If Pis to 
be even, then P’(0) must be zero, and P(Q) cannot then be zero. If P is to be 
odd, then P(0) must be zero, and P’(0) cannot be zero. If both P and P’ are 
zero anywhere, then, as you can see from the numerical procedure, P’”’, P’, and 
P will be zero everywhere. The ‘‘correct’’ non-zero values for P(Q) (in the 
even case) or P’(O) (in the odd case) would be given by the normalization 


condition f |P|?dx= f P2dx = 1. The easiest thing to do, however, 
OO 


is to make P(O) = 1 (or P'(Q) = 1), calculate the right shapes and energies for 
the states and, then, make the normalization correct. After all, as you can see 
from the Schrodinger equation, if P is a solution, then any constant times P is 
also a solution. We can normalize the wave function after we find it. This 
theorem has allowed us to set the initial values of P and P’ at x = O and then 
integrate out to large x before examining the behavior to the tail of the wave 
function. There are other ways to start the procedure: that is, there are other 
ways to initialize P and P’ at some known x. The procedure we’re using works 
for all symmetric potentials. 


THE HARMONIC OSCILLATOR 


Let’s work through a representative exercise. In fact, let’s start with a problem 
that analytic methods can solve—the harmonic oscillator. The potential for a 
harmonic oscillator is V(x) = kx2/2 where k is the so-called “spring constant.”’ 
The prototype of a harmonic oscillator is a mass oscillating on a spring, but 
we're actually treating the problem of an electron in a potential proportional 
to x?. For definiteness, let’s set the spring constant, k, equal to 1, so that 
V(x) = x?/2. 


The harmonic oscillator can be solved analytically. The result is that the wave 
functions are Hermite polynomials and that the allowed energies are E, = 
(n+1/2)ha where w = V/k/m. tn atomic units E, = (nt+1/2)w@ which means 
that, for our problem, E, = (n+1/2) or E = .5, 1.5, 2.5, . . . Each energy 
has a wave function that obeys the boundary conditions; no other energies 





will work. Since these energies are known, we can check the numerical cal- 
culations by seeing if these energies work and others don't. Furthermore, if 
some other potential, V(x), looks sort-of like a harmonic oscillator potential, 
then these energies can also be used as first guesses at the energies for 
that potential. 


SAMPLE SOLUTION 
1. Plot the potential. 


It’s always a good idea to get a feeling for the potential the particle 
experiences. You can do this by a short program which prints out 
values which you then plot on graph paper. One such program (with its 
RUN for this problem) would be: 


POTENT 


18 FOR X8e@ TO 3 STEP -2 
20 PRINT X@sXO*X6/2 


30 NEXT X86 

46 END 

END 

RUN 

POTENT 

@ 8 

22 82 
eA ¢ 68 
06 «I8 
08 32 
1 5 
1-2 72 
12a «98 
126 1-26 
18 1-62 
2 2 
2:2 2-42 
2e4 2-88 
2-6 3< 38 
2<8 3092 
Je 4-5 


END 


—_ 
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You can also use the terminal plotting subroutines discussed in the 
Appendix. The program changes you write (and a RUN for this prob- 
lem) could be: 


S§ READ X8sX9sY8,Y9,Z28 

6 DATA 8.328, 5,1 

7? GOSUB 9666 

19 FOR X@=@ TO 3 STEP -.2 
16 LET YOeX8exe/2 

26 GOSUB 391868 

38 WEXT X6 

35 GOSUB 9266 

49 ED 


APPEN D- TTY PLO 
RUN 
POTENT 


MAX Y= 5 
XXXXKXXKXKXAXKXX KXAAX XK AKA KX AK KK XX KX KAAKKK RAK KX MA AKKK XXKAAK KK AXA AK KAKA AKA KKK 


ROKER KKK KERR KEK KKK 


i i j 
AMXXKXKXAXAXXKAXX AKKXKXAKAAX KA KXKKXX KAKAK XX KAR XXKRK NKR KX KX KANK KARR KKK KK AK X AKAKA KKK XX 
MIN Y= @ MIN X= @ MAX X=" 3 


The potential is, of course, just a parabola. It is symmetric (V{—x) = 
V(x)} so we can use our computer approach as discussed above. 


2.(a) Sketch the wave function for the lowest energy (the ground state) 
of the potential V(x) = x2/2. Now find the wave function and energy of that 
state with the computer. 


Let’s use a program like that discussed above with the defined function 
for the potential. Then we just try different energies until the wave 
function behaves at large x. Large x usually means x several times the 
classical turning point (where E = V(x)); for this lowest state “large x’’ 
(the end of the interval of integration) can be 3 or 4. | 


Usually the lowest state for a symmetric potential is even; higher states 
usually alternate odd, even, odd, even, .. . A sensible expectation is to 
try for an even ground state. Notice, if the lowest state weren't even, 
the lowest odd state would just come out with a lower energy. The 
program and a RUN look like: 


tDSCHR 
1@ DEF FNvcx>=x«x72— Potential, V(x) 
20 PRINT “END OF INTERVAL, DELTA-X» PARITY (@@EVENS !=0DD>, ENERGY 7” 
38 INPUT X9eX7, PE 
46 LET PG=P1=X69¢ 
S@ LET X6"x9/20 € es ae 
66 IF Pel THEN 96 Initialization 
78 LET Peel 
88 GOTO 186 
9@ LET Pi=1 ; 
199 LET PReleCFNVCXEOXT /2)~£)e( POSP 16X72) : 
118 LET P@@Pee( PleP20X7/2) 4X7 Integration of the 
126 LET PlePl+P2exT Schrodinger equation 
136 LET X@eX@eX7 
14@ IF X@<X6 THEN 178 : 
1S@ LET X60X6+X9/28 Print loop 
168 PRINT XO,P6, POePa 
176 IF x@<x9 THEN 106— Aeturn for next Ax 
18@ PRINT 
196 Goto 2g— Aeturn for new parameters 
268 END 
RUN 
1DSCHR 
END OF INTERVAL, DELTA*%» PARITY (8#EVEN3 1*0DD), ENERGY? 
730 -e81s 0s 2&5 
01S ©9899 34 «9 799 69 
031 095782 0917419 
+ 46 0989852 082783 
261 2847857 °718861 
076 0775826 © 6819986 
2989999 + 698172 ° 487448 
1006 © 619367 ° 383615 
1e2! «543639 © 295544 
1236 ¢ &7 4765 - 2254862 
125! 2 418995 173652 
1-66 0378122 013699 
te8} 339714 ° 135406 
1096 © 327534 2187279 
2-1 © 33576 2112738 
2.25 © 379875 © 137326 
2. 48001 44853 © 19.4067 
2.5506! «559 386 0312913 
2-7¢6661 749964 2562476 
2085662 1.5051 1619356 
3080002 1-52538 2.32678 
END OF INTERVAL, DELTA-X» PARITY (@@EVENJ 1®0DD). ENERGY? 
72322081, 8,-55 
o!S 2987692 0975536 
231 094836 0899386 
046 2889394 791821 
061 0812719 © 668512 
°76 °722762 » 522385 
2989999 + 624298 + 389748 
1-86 «5219 66 e 272449 
1221 © 419833 © 17626 
1636 032104 « 183866 
1051 © 227348 5.1778 1E-@2 
1066 2139977 1-95935E-02 
1081 5.74888 E-682 3- 3949 TE-63 
1096 ~ 20 23332E-62 A.987TIE-84 
21 -9 8 128 6E-@2 9 620 22E-93 
2.25 ~» 186798 3.489 35E-82 
2. AGGG! -. 292586 8 5559 6E-@2 
2. $5621 ~~ 429581 018 4a7t 
2078801 ~. 619735 - 384071 
2.85002 ~.898843 088648 1 
3. eea¢2 -1032131 1.74585 
END OF INTERVAL, DELTA=X» PARITY (@#EVEN} i=0DD), ENERGY? 
23e0e8ls Gee 
015 986813 0977753 
°31 2953086 0988373 
246 -899 685 - 889238 
61 830232 » 689 285 
076 0749 162 561244 
°989999 © 6689 68 + 436878 
1-86 0578181 - 325167 
1.2! © 486921 0231265 
136 ° 39 6689 + 157299 
1eS1 2319881 + 162272 
1266 0252129 6¢ 3869 1E-82 
181 © 194383 3077732E-82 
1396 © 146482 2821487 
Ze1 211624 le 21528E-82 
2025 © 879544 6¢32724E-83 
2. 40001 5. 611 28E-82 3- 148 6SE-93 
2.55001 308 OV ATE-B2 16 49728E-63 
2078061 2e 6077 1E-82 6e8G016E-64 
2085082 1671633E-82 2094579 E-684 
3080602 Le 10143E-62 1021314E-64 
END OF INTERVAL, DELTA-X. PARITY (@*EVENS 12#0DD), ENERGY? 


? 
END 
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As we expected, the energy of the state is near .5. You can get higher 
accuracy by using a smaller step size (Ax), but because the computer 
keeps numbers to only finite accuracy and because the calculational 
method is not perfectly accurate, the computer may not give you ex- 
actly .5 as an answer. The difference from .5 is a measure of the accuracy 
of the method. 


2.(b) In which regions of x is the probability of finding the particle nearly 
zero? In which regions of x is the probability of finding the particle relatively 
large? 


Since P is real, | P12 =P2. By looking at the printout you can see there 
is one loop (or anti-node) in P? and there are no nodes. The probability 
is near zero only for large positive and large negative x. The probability 
is Jarge near x = 0. You can also plot the wave function and its square 
on the terminal. The changes you make and the RUN are: 


5 LET X8#Y8=6 
6 LET Y9u! 

35 GOSUB 9668 
168 LET Z6s1 
162 LET Y@=Pa 
163 GOSUB 9188 


164 LET Z6@=2 

165 LET Y@=ePOePe 
166 GOSUB 9188 
188 GOSUB 9288 
APPEN D- TTY PLO 
RUN 

IDSCAR 


END OF INTERVAL.» DEL TA-X» PARITY (O=EVEN3 180DD). ENERGY? 

73Jo281,82-5 

Max Ye 1 

KOMIK AA HAKAN AIA AKAN KKKKAKK AK MAK AMMMAM AMMA MAAK MARMARA AMA MAMMA MM A MMMM 


2 

2 2 : 1 i 
2 2 2 2 2 2 
XAXKXKKAXAAXAAKAKAAKAAKAAKRAAKAA AKA AKA KAKMK AKA A HAMANN HEAR NHN MMH HA MAAK AAA KAR AKAN 
MIN Y= 6 MIN Xs 6 MAX X® 3 


END OF INTERVAL, DELTA=Xs PARITY (@=@EVEN3 1#0DD), ENERGY? 
? 


GE me ee eee ee RK 


END 





2.(c) Find the normalization constant for the wave function of this state. 


Coo 
To normalize a wave function, you calculate A = f |P|7dx. Since 
— oO 


P(x) is real, | P|? = P*. The normalization constant (i.e., the number 


. ; oo 
you multiply P by everywhere) is then T/JA. The total integral f 


—_0O 
P?dx is just twice the sum of P2Ax from x = Oto your end of interval. 


Thus, with a minor change in your wave function program, you can 
calculate A directly. 


45 LET P32 

115 LET P38P3+ Poe Ppoex7 
168 PRINT X8,P@,P0«PO, P3 
RUN 

1 DSCHR 


END OF INTERVAL, DELTA+*Xs PARITY (8s EVENS 120DD), ENERGY? 
230+ Bl,Bs-5 


°15 988813 977751 ° 148 769 
e3t °953886 998373 ° 299887 
246 » 899 605 «809268 » 428555 
- 61 «838232 ° 689265 ° 548528 
#76 749162 »$61244 + 633699 
989999 e 6689 68 ¢ 436878 «787828 
1+@6 ° 570181 325107 764219 
1221 « 488921 °231285 °88S237 
1-36 39 6689 ¢ 157299 °833765 
1251 3196881 © 182272 °8527358 
1-66 252129 62-3569 1E-82 «864797 
1.81 © 194353 3¢77732E~ 82 °87213 

1696 «146482 © 821457 876393 
2-1 ellgea 1e21528E-82 878645 
2-25 © 879544 66 32724E-03 87996 

2¢ 48681 §$-61128E-82 3e 148 6SE-83 °88863 

2-55861 3-869 47 E-82 f- 49728 E-83 -888957 
2-78661 2- 68077 1E-82 6-88B16E-84 8811689 
2.858¢2 1-7 1633E<B2 2-9 AST9E-84 -881176 
3+ 88862 ls 181 43E-92 le21314E-64 +88 1204 


Pe The normalization constant is 1/4/2(.881) = .753. 


2.(d) Find the position x, such that the probability of finding the particle 
between ~x, and +x, is 1/2. 
x 


0 
What you really want is the x, such that f | P| 4dx = 1/2 whenP is 


~X, 
normalized. All you need to do is look down the printout from Part 2(c) 
to find the point where the sum of P2Ax = A/2. The point is very near 
x = .47 atomic units = .150 A. (For higher accuracy you could change 
the print-test so values are printed more often on the terminal. You 
might also plot the P2Ax function and estimate the point x, graphically.) 


3.(a) Sketch the wave function for the second lowest energy of V{x) = x?/2. 
Now find the energy and wave function for this state with the computer. 


3.(b) In which regions of x is the probability of finding the particle nearly zero? 
In which regions of x is the probability of finding the particle relatively large? 


3.(c) Find the normalization constant for this state. 
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3.(d) Find the position such that the probability of finding the particle be- 


tween “Xo and #X is 1/2. 


The RUNs for this part of the problem follow: 


i 


END OF INTERVAL, 
PAegeBle ls 1-48 


END OF INTERVAL. 
2A, e6le121-52 


END OF INTERVAL, 


? 


END OF INTERVAL» 


? 


DSCHR 


2! 
24} 
«61 
©869999 
ie81 
1-21 
be 4i 
1-6! 
te81 
2061 
202 
2. 49081 
2+ 6808! 
2088882 
3+« 86862 
3¢ 26883 
3e 40003 
3+ 66663 
386664 
4. G6604 


2! 

241 

« 61 
2869999 
1-61 
1e21 
104] 

le 61 
1281 

2e 81 
2.2 

2- aG861 
2- 60601 
2886662 
3- 86262 
3e 28863 
3e 46883 
3e 68663 
3086004 
4-80604 


A, eo Ble la le5 
2! 

24! 

e 61 

» §69999 
Te@l 
1e2i 
1241 
1261 
1e81 
2081 
2e2 

2- 4886) 
fe 68901 
20690082 
3e BEBS2 
3e 26883 
3<¢ 48003 
3« 66003 
3¢8 8884 
4e F8684 


END 


¢ 283415 

e 37 69 38 

e 566427 

e 58 3447 

¢ 66645 
581894 
521783 

« 44949 
351766 

° 2666 

«1956 
13469 1 

8 e848 11E-22 
e 85549 6 
«833239 
816988 
1e«@2688E-82 
§- 1 LISSE-83 
2e O@2241E-03 
“1e31LIQ@E-84 


DEL TA-X~s 


DELTA=X, PARITY (@*EVEN3 
«285476 As 22204E-62 
° 377375 + 142412 
° 587786 ° 2578 46 
« 58 6386 23438 48 
© 618624 0374084 
° 58987 ©3479 47 
253362 -28alti 
© 455398 «287388 
¢ 37085 ° 13753 
2298693 8 + ASO2LE-92 
© 225635 5. 1O81SE-98 
+ 174964 3. O2774E-82 
o 141619 «626656 
©139921 10 71403E-82 
© 146246 2. 13878E-62 
2197919 3.927 18E-62 
© 309 367 9257879 E°62 
0531873 + 28 2889 
°977555 0955614 
[+89 648 3059 66S 

DELTA-X» PARITY (@=EVEN} 
2285354 4e21783E-82 
- 37658! ~ 141753 
«505872 © 255897 
© 588517 0337 
© 681299 - 361561 
° 57397 0329441 
° 5106S © 266763 
+ 425769 «181279 
. 332998 - 118888 
© 243817 $.908874E-62 
«166152 20 76Q65E~82 
96609 1E-82  9233332E-63 
307 26 68E-2 1639631E-93 

-16 67988E-02 2.82199E-84 

-7. 4578 6E~82 $.56196E-83 

-2 151841 2.28 133E-82 

= 272922 7+ 48 GTE-82 

-. 492861 ~242124 

=0916188 0839253 

91678854 3.17833 


4e219S53E-82 
« 142682 
« 256469 
° 34941! 
+ 367782 
¢ 338 691 
© 272257 
194831 
«123739 
7e18753E-62 
308 259 5E-8e 
1-8 1416E-82 
7-8 289 8E-83 
3° #798 LE-23 
1» 1848 3E-83 
30 6039 3E-24 
le BS429 E-84 
2e 613 18E-05 
&- 89 812E-86 
le 7T26S54E~- 68 


PARITY (@= EVEN) 


1s ODD). 


12 0DD).» 


DELTA-X%s PARITY (@#EVEN) t#ODD), 


l* ODD). 


ENERGY ? 


ENERGY ? 


EN ERGY 7 


ENERGY? 





3.fe) 


3.(f) 
3.(g) 
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To plot the functions, start again with the basic program and add: 


§ READ X8sY8sY9 
6 DATA 82 Bs 1 

35 GOSUB 9882 
168 LET Z9#1 
1631 LET Y@=ePr@ 
162 GOSUB 9162 
163 LET Ze@e#2 
164 LET YOaPS*Pg 
165 GOSUB 9182 
188 GOSUB 92862 
APPEN D- TTY PLO 
RUN 

IDSCHR 


END OF INTERVAL» DELTA=X%s PARITY (8#EVENS {!=0DD),s ENERGY? 

240 eBte loa leS 

MAX Y® 1 

XXXKKKK XKXK KXKXKKK KKK KK KX KKXAKAX KX RX KKAXKAKKKRNRAK KX KX KKAKKXAKK KKH AK K ARXK KX KK AK 


Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 11 

Y i 

Y 1 

Y 1 

¥ 

Y ! 

Y 

Y 1 2 1 

Y 2 2 

Y 

Y 2 2 1 

x 

Y 1 2 1 

Y 

Y 2 2 i 

Y ! 

Y 2 1 
Y e 2 2 11 
Y @e2 22 22 e2 
XAXAXKXXKAXAAAAAKAAXAKARAKAAKK KAAAAAAXNAANK AANA ANA AXA HAHK KXAN K ANKXAX XK X 
MIN Ys @ MIN X= @ MAX X= a 
How many nodes are there in [P{x)]?? 


There are 2 loops and 1 node between —co and too, The rule of thumb 
is that the lowest energy state has O nodes, the next state has 1 node, 
etc. The probability is nearly zero around x = O and for large positive 
and negative x. The probability is large near x = +1 and x = -1. 


Find the normalization constant. 


Find the point xX, such that the probability of finding the particle be- 


tween -X, and Xz, is 1/2. 


Starting again from the basic program, add: 


45 LET P39¢ 

115 LET P3eP3+peepaex7 — adds up P2 Ax 
160 PRINT X@,P@,P@#P@, P3 

RUN 

1DSCHR 


END OF INTERVAL, DELTA-Xs PARITY (QO=EVEN] 12 0DD), ENERGY? 
PAs e- Bl, 4,165 


21 » 2865415 4-21953E-62 3« 22862E-83 
«4} » 37 6938 e 142682 2-15839E-82 
«61 © 586427 «256469 6¢ 20956E-82 
°889999 © 58 3447 ° 348411 123885 
1o@1 « 68645 ° 367782 «194978 

e211 «581894 33868! ° 266289 
tea) 521783 2272257 «327448 

Le 61 «44849 194631 © 37369 

1-8! 351766 + 123739 «494872 
2-01 « 2666 7e!87TS53E-62 © 423768 

2.2 ©1986 3e8259SE-82 © 433746 

2e 489001 © 13469 1 1¢8 1416E-82 © 439 88 1 

2- 68681 S-84611E-62 728289 6E-83 2441583 
2-880682 * 8$5496 393-8796 1E-83 © 442595 

3+ 88682 833239 le }8@83E-83 © 442883 

3+ 28883 «818984 3-¢ 6839 3E-84 © 443914 

3+ 46883 te 82688 E-82 te 85449 E- 64 © 443854 

3+ 66883 Se 11185E-83 2e 61IIBE-B5 » 443866 
3-88004 2- 62241E-83 4-896812E-96 «443868 
4288064 “1-31169E-94 1-72054E- 88 ° 443968 


Note: Curve 7 is wave function x = 0 tox =4; 


curve 2 is the square of the wave function. 
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The normalization constant is 1/\/2(.443) = 1.062. The point at which 
the sum of P2Ax = .221 is approximately x = 1.1 atomic unit = .532 A. 


4, Find the difference in energy between these states in eV, Hertz, and A. 


5. fa) 


5. (b) 


E = 1 atomic unit = 27.210: eV = 6.58 x 1015 Hertz; lambda = 456 A. 
(Using 1 eV = 1.602 x 10°12 ergs; h = 6.624 x 10°?” erg-sec; c = 
2.998 x 1019 cm/sec) 


Find the behavior of each state near x = O. 


Consider 2n(P(O)-P) versus &n(x) for small x. What is the power de- 
pendence of P for small x? 


The slope of a straight line on a log-log plot is the power, B, to which x 
is raised in (P(O)-P) ~ x®. You may want to plot your results on graph 
paper although you can calculate the result directly from printout. 


Starting again from the basic program: 


GET~ 1DSCHR 

148 

1$8 

160 PRINT LOG(X@).LO0GC 1-P9) 


17@ IF X@<el THEN 188 

RUN 

1DSCHR 

END OF INTERVAL.» DELTA~X» PARITY (O=sEVEN} 1t< ODD), ENERGY? 
730 0¢O15 0005 

~4. 69517 -9.98451 

~3091282 -8-51762 

-3- 58656 -7 «78676 

~3e 21888 -7213148 

-2099573 - 6« 68535 

2661341 - 6032893 

~2- 65926 ~6e91294 

~2652573 #5074623 

-2- 48795 5.51187 

-2-39259 -5- 38883 

END OF INTERVAL, DELTA-X» PARITY (@#EVENS i2O0DD). ENERGY? 


? 
END 


So this wave function comes off of its value at x = O (P(O)=1) quad- 
ratically. The normalization doesn’t change this fact. The first exicted 
state (the second one found above) comes out of P(Q) = O linearly in x. 
Any even state must come off its x = 0 value as an even power of x; any 
odd state must start out from x = O as an odd power of x. 


You actually know a good deal about a wave function before you start 
a problem. Often the Schrodinger equation can be approximated in 
some regions of space by an equation with known solutions. The idea 
of these asymptotic solutions for small and large x will be important 
again in Chapter two. 


Find the behavior of each state for large x. 


Consider &n(P) versus x for large x. Is the dependence of the wave 
function beyond the classical turning point faster than e%* where a is 
some constant? ; 


'f the potential became a constant, VO, for large enough x, the tail of 
the wave function would approach ex p(-ax) where a = +/2(VO-E). Many 
notentials do approach constant values, but V(x) = x?/2 clearly ap- 
proaches infinity. So, you expect the wave function to tail off to zero 
faster than exponentially. You may want to plot the results on graph 


paper. 
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Starting again from the basic program, change: 


140 IF X@<2-e9 THEN 170 
168 PRINT X@sLOGC PO) 


RUN 

iDSCHR 

END OF INTERVAL» DELTA-Xse PARITY (@@EVEN3 JI=#O0DD)s ENERGY? 

230e8tsBeeS 

2098882 ~4- 2699 
2091862 #4223922 
2092002 “Ae 268 67 
293882 ~4- 29823 
2054982 ~4- 32791 
2095862 “4.35771 
2096882 - 4.38763 
2.97882 ~4. 41768 
2098882 ~4. 44735 
2.99662 4-473 14 


So the tail of the wave function falls off faster than exponentially be- 
cause V(x) is still rising fast. The normalization doesn’t affect this re- 
sult either. Interestingly, if you assume that the tail behaves as 
exp({-(a@(x)x) where a(x) = /2(V(x)-E}), you can get good agreement 
with the asymptotic behavior of P at large x. 

EXERCISES — Use computer when appropriate 

1. For the potential V(x) =x?: 
1. Plot the potential forO<=x <= 85. 
2.(a) Find the ground state energy and wave function. 
2.(b) How many nodes are there in [P(x)] 2? 


2.(c} Find the normalization constant. 


2.(d) Find the point Xo such that the probability of finding the particle 
between -X, and xq is 1/2. 


3.fa) Find the first excited state wave function and energy. 
3.{b) How many nodes are there in [P(x)] 2? 
3.fe) Find the normalization constant. 


3.(d) Find the point Xg such that the probability of finding the particle 
between -x,and x, is 1/2. 


4, Calculate the energy difference between these two lowest states, AE, 
in atomic units, eV, Hertz, and A. 


5.(a) Find the behavior of each state near x = 0. 

5.(b) Find the behavior of each state for large x. 

Hint: The energies are known from the harmonic oscillator with k = 2. 
Il. For the potential V(x) = (x2/2) + 4 exp(-x?), which is a harmonic oscil- 
lator with a Gaussian bump in the center (a model of the potential of the 
nitrogen atom in an ammonia molecule): 


1. Plot the potential forO< =x <= 85, 


2.(a) Find the ground state energy and wave function. 


14 UO QUANTUM MECHANICS 


2. (b) How many nodes are there in [P(x)] 2? 
2.(c}) Find the normalization constant. 


2.(d) Find the point x, such that the probability of finding the particle 
between ~x, and x, is 1/2. 


3.{a) Find the first excited state wave function and energy. 
3.(b) How many nodes are there in [P(x)] 2? 
3.(e) Find the normalization constant. 


3.(d) Find the point xq such that the probability of finding the particle 
between -x, and X, is 1/2. 


4, Calculate the energy difference, AE, between these two lowest states 
in eV, Hertz, and A, 


5.(a) Find the behavior of each state near x = 0. 
5.(b) Find the behavior of each state for large x. 
6. Compare the energies and wave functions of these two states to those 
discussed in the chapter {i.e., the two lowest energy states for the har- 


monic oscillator x2/2). 


Hint: Use the proper harmonic oscillator energies as first guesses and 
then hunt. 


itl. For the potential V(x) = 16x?/(x2+16) which starts out as a harmonic 
oscillator but then becomes constant for x >>4" 


1. Plot the potential forO<=x <= 5. 

2. (a) Find the ground state energy and wave function. 
2.(b) How many nodes are there in [P(x)] 2? 

2.(c) Find the normalization constant. 


2.(d) Find the point x, such that the probability of finding the particle 
between -x, and X4 is 1/2. 


3.(a) Find the first excited state wave function and energy. 
3.(b) How many noises are there in [P{x)] #? 
3.{c) Find the normalization constant. 


3.(d) Find the point xg such that the probability of finding the particle 
between -X, and Xq is 1/2. 


4. Calculate the energy difference, AE, between these two lowest states 
in eV, Hertz, and A. 


5.{a) Find the behavior of each state near x = 0. 
5.(b}) Find the behavior of each state for large x. 


6. Compare the energies and wave functions of these two states to those 
discussed in the chapter {i.e., the two lowest energy states for the har- 
monic oscillator x2/2.) 


Hint: Use the proper harmonic oscillator energies as first guesses and 
then hunt. 
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IV. For the potential V(x) = - 10 exp(-x?/2), which is a Gaussian well: 


1. Plot the potential forO<=x=<=85. 

2.(a} Find the ground state energy and wave function. 
2.(b}) How many nodes are there in [P(x}] 2? 

2.{c) Find the normalization constant. 


2.(d) Find the point x, such that the probability of finding the particle 
between -x, and Xq is 1/2. 


3.(a) Find the first excited state wave function and energy. 
3.({b) How many nodes are there in [P{x}] 2? 
3.{c) Find the normalization constant. 


3.(d) Find the point x, such that the probability of finding the particle 
between -x, and xq, is 1/2. 


4, Calculate the energy difference, AE, between these two lowest states 
in eV, Hertz, and A. 


5.(a) Find the behavior of each state near x = 0. 
5.(b) Find the behavior of each state for large x. 
6. Compare the energies and wave functions of these two states to those 
discussed in the chapter (i.e., the two lowest energy states for the har- 


monic oscillator x2 /2), 


Hint: Use the proper harmonic oscillator energies as first guesses and 
then hunt. 


*V. For the potential V(x) =-exp(-|x|/10)/|x|, which is the one-dimensional 
equivalent of the screened Coulomb potential (important in many electron 
atoms): 


1. Plot the potential forO< x <= 5, 

2.(a} Find the ground state energy and wave function. 
2.(b) How many nodes are there in [P{x)] 2? 

2.(c) Find the normalization constant. 


2.(d) Find the point x, such that the probability of finding the particle 
between -x, and Xx, is 1/2. 


3.fa} Find the first excited state wave function and energy. 
3.(b} How many nodes are there in [P(x)] 2? 
3.(c) Find the normalization constant. 


3.(d) Find the point x, such that the probability of finding the particle 
between -X,_ and x, is 1/2. 


4. Calculate the energy difference, AE, between these two lowest states 
in eV, Hertz, and A. 
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5.(a) Find the behavior of each state near x = 0. 


5.(b) Find the behavior of each state for large x. 


Hint: Because the potential is infinite at x = 0, only odd states are 
allowed. The (unscreened) Coulomb potential (the hydrogen atom) 
energies for these states are -.5 and -.125; you will find both energies 
moved up. This potential is important in three dimensional quantum 
mechanics (Chapter Two). There, it turns out, the wave function is 
Rir) = P(r)/r where P(r) is the wave function you derive in this problem. 
You might look at the behavior of P(r)/r in your printout. 


For the potential V(x) = (|xl-1)2/2, the biharmonic oscillator: 
1. Plot the potential forO<=x <= 85. 
2 (a) Find the ground state energy and wave function. 
2.(b) How many nodes are there in [P(x)] 2? 
2.(c} Find the normalization constant. 


2.(d) Find the point x. such that the probability of finding the particle 
between -X, and Xz fs 1/2. 


3.(a) Find the 4th excited state wave function and energy. 
3.(o) How many nodes are there in [P{x)] #? 
3.({c) Find the normalization constant. 


3.(d) Find the point x, such that the probability of finding the particle 
between -x, and Xq is 1/2. 


4, Calculate the energy difference, AE, between these two lowest states 
in eV, Hertz, and A, 


5.(a) Find the behavior of each state near x = 0. 


5.({b) Find the behavior of each state for large x. 


6. Compare the energies and wave functions of these two states to those 
discussed in the chapter {i.e., the two lowest energy states for the har- 
monic oscillator x?/2). 


Hint:Use the proper harmonic oscillator energies as first guesses and 
then hunt. 


For the potential V(x) = x?/2: 
1. Plot the potential forO<=x <= 5. 


2a) Find the 3rd excited state (the fourth state) energy and wave 
function. 


2.(b) How many nodes are there in [P(x)]?? 
2.(c}) Find the normalization constant. 


2.(d) Find the point x, such that the probability of finding the particle 
between -x, and Xq is 1/2. 


3.{a}) Find the 4th excited state wave function and energy. 
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3.{b) How many nodes are there in [P{x)] 2? 
3.(c) Find the normalization constant. 


3.(d) Find the point xq such that the probability of finding the particle 
between -x, and xX, is 1/2. 


4, Calculate the energy difference, \E, between these two lowest states 
in eV, Hertz, and A. 


5.(a) Find the behavior of each state near x = 0. 
5.(b) Find the behavior of each state for large x. 
6. Compare the energies and wave functions of these two states to those 
discussed in the chapter {(i.e., the two lowest energy states for the har- 


monic oscillator x2/2). 


Hint: Use the proper harmonic oscillator energies as first guesses and 
then hunt, 


**VIIl. The family of potentials, V(x) =| x |™ for different integers, m. form 
an interesting progression, the m = 2 member is a harmonic oscillator with 


energies, E, = (n + 1/2)./2; the limit as m > © is the infinite square well. 


1, Find the wave functions and energies for the ground state with m=2 
4, 6, 8, and 10. 


2. Find the wave functions and energies for the first excited state with 
m= 2, 4,6, 8 and 10. 


3. Plot the energies as functions of the power, m. 
4. Interpret the shapes of the curves from Part 3. 
"*IX. For the potential V(x) = (-x2/2) + A exp(-x2): 


1, Find the energies and wave functions of the ground and first excited 
states as functions of the height, A, of the Gaussian bump. 


2. Plot the energies and their difference versus the height of the bump. 


3. Explain the changes in shapes of the wave functions. What conclu- 
sions can you draw about the energies? Consider the extrapolations to 
A = 0 and A = -% consider the relationship of the kinetic and potential 
energy terms in the Schrodinger equation. 
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CHAPTER TWO: THREE-DIMENSIONAL BOUND STATES 
INTRODUCTION 

In genera!, three-dimensional quantum mechanics is thought to be much 

more difficult than the one-dimensional case. For spherically symmetric po- 

tentials, the computer makes three-dimensional quantum mechanics just as 


easy as One-dimensional. 


The three-dimensional time-independent Schrodinger equation is 


h2 
5m V2 WIT) + VIF) WP) = EV (r) (4) 


where W is the wave function for the energy E and V(r) is the potential. 
In polar coordinates (r, 6, d) this equation becomes 


W243 2 ay) aad (sino 28) +, FY] 
~ 2m r2 at dr_| mr? [sind a0 (sino a0) * sin26 0p? (5) 


+V(r,0,¢@)¥ = Ey 








lf V(r, 8, o) is spherically symmetric, i.e., Vir) = Vir), then Equation (5) 
separates in polar coordinates. The wave function can be written as 
Wy = R(r)Y(@, ¢). You substitute this product in for W, put all the terms 
containing R(r) on one side and all those containing Y(@, $) on the other; 
then set both sides equal to a constant (&(2+1) where 2 is any integer that is 
convenient). This provides the following equations for R(r) and Y(@, ¢): 


2 2 
#4 4[ ash] + vin + We Rir) = ER(r) (6) 
1 Y (0 1 392Y(0 
sind > sin oe Teneo ae + 2(2+1) Y (2,4) =0 (7) 


THE EQUATIONS IN @ and ¢ 


The equation for Y can be separated again as Y(@, ¢) = O(8)d(¢) with the 


result that 
1 df., do m? a 
aa as oo) * (eee : a) =a i 
d2@ 
doz m2p = 0 (9) 


The ®(¢) equation has the solutions 


@(d) = eiM® and etime (10) 


where m is an integer such that | mj < = &. Since we will only use ® in 
applications involving | y |2 for the ¢ = 0 -¢ = 180° plane, we don’t have 
to worry about ® at all. If you need to use ®, you can call the SIN{ ) and 
COS() functions on the computer. 
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The equations for ©(0) is the associated Legendre equation; its solutions are 
known as associated Legendre polynomials, £5", so 


O(0) = £6" (cos A) (11) 


The other solution to this second order differential equation for @ isn’t finite 
throughout the range of @ and so is not allowed to occur in most physics 
problems. So the solution for the angular parts of the full wave function, y, is 


Y (6,0) = No, £5" (cos) gt ime (12) 
where No, is just a normalization constant. 


You can calculate values for ©(@) and P(d) easily on the computer. To cal- 
culate © you could integrate the differential equation numerically. However, 
like sines and cosines, associated Legendre polynomials are well Known and 
can be looked up in books. The first few angular wave functions, Yg" (@, 6) are 


se = 1/741 


Ys =¥ = cosé , v = = sind et iM? Y;' = Jf = sind eM? 


One possible program to calculate the Yo. s is the following: 


16 PRINT “"L (@ OR 105 M (C/M/<8L) 2" 
26 INPUT LM 
38 LET Pil#3-14159 
46 LET N@@=1/SORCaePI> | 
SO LET Be | 
6@ LET N2=SOR(3/(8*P1)) 
76 FOR Tl=@ TO Pl STEP P1/6~<—Step across 0 

88 LET CleCcos(TID 

9@ LET SimSINCTI) 

1@@ PRINT “THETA ="3 160¢TI/P1 

116 FOR T2e@ TO @ePl STEP 2«Pt/6 «—Step across 
186 PRINT 160«T2/P1, 

136 IF L>@ THEN 166 0 

146 PRINT N@ Y, 


Useful constants for calculation 


186 GOTO 236 0 
166 IF @ <> © THEN 196 | 0 
17@ PRIBT N18Ci | 

189 GOTO 23¢ ( “7 


190 IF Mao] THEN 2e8 Pang 
208 PRINT N2eSieCOS(T2>oN2eSiaSIN(T2)> Y 

218 GOTO e38 , 7 

220 PRINT N2¥SI*COS(T2),-2eS1¥SIN(T2)—@__y, -7 
238 NEXT T2~< 6 1 


278 gore 18<— Return for new parameters 
D 


THE RADIAL EQUATION 


The real problem of three-dimensional, spherically-symmetric-potential, quan- 
tum mechanics is to solve the radial wave equation, equation (6). Only one 
problem (the hydrogen atom, which is the unscreened Coulomb potential) 
can be solved analytically. All other problems need approximate methods of 
one sort or another (perturbation theory or variational methods), and sooner 
or later even these approximate methods need a computer. 


Using the computer, an entirely different approach can be taken. You can 
integrate iteratively to find the radial wave function and the energy. Such a 
procedure can be carried out for any potential, Vir). Just as was true in 
Chapter One on one-dimensional bound states, the computer approach for 
three-dimensional bound states is easy to understand but very general. 
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First, we reduce the three-dimensional radial wave equation to one that looks 
just like the one-dimensional Schrodinger equation discussed in Chapter One. 
Let P(r) = rR(r). You then plug P(r) into the radial wave equation (6) and 
find the equation P(r) satisfies. The result is 


2 d2P(r) h22(2+1) 


In atomic units (for which h=m=e=1; see the Preface) this equation becomes 








1 d2P(r) Q(2+1) 
oa vee ooo ]Pc = EP(r) (14) 
or 
" Q(L+ 1) 
P = 2] vir + Se |r (15) 


This last equation looks just like the one-dimensional Schrodinger equation 
discussed in Chapter One. The only differences are (1) that you are integrating 
with respect to r instead of x (so, for example, the range of r is 0 to +°°) and 
(2) that the potential V(x) is replaced by the effective potential V...(r) = V(r) 

Q(L+1) 
+ 

2r2 

Schrodinger equation. 





. We can solve this equation just as we did the one-dimensional 


Suppose you know the reduced radial wave function, P(r), and its first deriva- 
tive, P(r) = a at some point r,. Suppose you also know the behavior of P 
- somewhere else (for example, (1) that P(r) must go smoothly to zero as r 
approaches zero (because R{r) must remain bounded)). Then, for any poten- 
tial, V(r), you can choose an energy E, set P(r) and P(r) at that first point fee 
integrate from r, to the second region, and examine P(r). If the energy you 
chose was an eigenvalue (one of the allowed energies for the system), then 
P(r) (and R(r)) will behave correctly. If E was not an eigenvalue, P(r) will 
diverge. The divergence will be to opposite infinities on opposite sides of 
an eigenvalue. 


Using a computer, the integration can be performed in a step-wise way. 
Knowing P(r) and P’(r) at some point r, you find the values at r+Ar from 


P(rtAr) = P(r) + PS) Ar (16) 


P'(rt+Ar) = P’(r) + P(e) Ar (17) 


d?P/dr2 can be calculated from the radial wave equation for P. The derivatives 
are evaluated in the middle of each interval [r, r+Ar] so that the calculations 
are more accurate (see the Appendix). 


Figure 2 shows a simplified flow chart for the strategy of the computer calcu- 
lation. Before we can discuss an implementation of this strategy, we must 
discuss the initialization of P and P’ further. 





Figure 2. Flow Chart of Computer Approach Three- 
Dimensional Quantum Mechanics 
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INITIALIZATION 


There are a number of ways to initialize P(r) and P’(r). Two ways are most 
common: the first initializes the (reduced) wave function, P, and its derivative 
at large r (and then observes the behavior at small r); the second initializes the 
(reduced) wave function, P, and its derivative at small r (and then observes the 
tail at large r). Large r in these contexts means several times the classical turn- 
ing radius (where the effective potential equals the chosen energy). Both 
methods involve asymptotic expansions of the wave function for some values 
of r. We'll discuss each method in turn. 


The first method expands the (reduced) wave function, P, and its derivative, 
P’ for large r. 1f the potential is such that, beyond some point, rg, the varia- 
Q(L+1) 
2r? 
say Vj. The Schrodinger equation for P takes on its asymptotic form 





tion of the effective potential is negligible, then V(r) + is aconstant, 


p” = 2(V, -E)P (18) 


whose solutions are exp(-ar) and exp(t+ar) where a=+/2(V,-E). Because there 
must not be infinite probability of finding the particle at infinity, the solution 
et is not allowed. The result is that, beyond the point rp, P(r)*exp(-ar) and 
therefore P’(r) © —aexp(-ar). You can initialize P and P’ far enough out in r and 
then integrate step-by-step back towards r = 0. The boundary condition 
imposed then at r = O is that R(r) must be bounded (hence P(r) must go to 
zero smoothly). Again the reason for this boundary condition is that infinite 
probabilities are not allowed physically. When you use this strategy, you will 
have to watch the behavior of P(r) and R(r) = P(r)/r carefully very near the 
origin (r=0). An implementation of the program using this initialization 
might look like: 


1@ DEF FNVCP)wageecre¢-12)-R1¢-6))— Potential, V(r) 

59 «PRINT “RCSTART)» RCEND) » DELTA-Rs ORBITAL #(L).» ENERGY?" 
3@ INPUT R9,R8,PR7+L,E 

a@ LET R7=-ABS(R7) 

5@ LET PO=EXP(~SQR( 26 ( FNVCR9)-E) ) ® POD sansa ase 

6@ LET Ple-SQR(2e¢FNVCR9)-E))*Pa Initialization at large r 
70 LET Reno 

8@ PRINT Re PO, P@/R 
98 LET RS=(RB-R9)/28 
1@6 LET R6=ReR7T/2 


118 LET Pemex  FNVCR6)4L4«(L+1) /C 2a R6*#R6)-E) *( POOP1 R772) Integration of 


12@ LET PO#Pe( P1+P2"R7/2)4R7 the Schrodinger 
136 LET P1=P1+P2#R7 

148 LET R=ReR? equation 

15@ IF R>RS THEN 188 { 


162 LET RS@RS+(RB-R9) /26 | Print loop 
L706 PRINT Re PO, PO/R | 

i8@ IF R>R8 THEN 10¢—Aeturn for next Ar 
198 PRINT 

20e GoTo 28— Aeturn for new parameters 
21@ END 
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The second method of initialization uses expansions of the reduced wave 
function, P(r), near the origin. Sometimes the Schrodinger equation reduces 
to a simple form for small r; sometimes the solution near the origin reduces to 
a solution already known. In either of these cases this second method of 
Initialization is the easier to use. A common example is the screened Coulomb 
potential, -exp(-r/r,)/r. This potential is a good approximation to the poten- 
tial experienced by the electrons in many electron atoms. For small r, this 
potential reduces to the unscreened Coulomb potential, -1/r. The analytic 
solution for the hydrogen atom (the unscreened Coulomb potential) can be 
expanded for small r, as R(r)r®. So the small r behavior of the wave function 
for any potential which reduces to -1/r for small r (in particular for the 
screened Coulomb potential) is given by: 


P(r) ~ r&t1 (19) 
so that 
P’(r) ~ (R+1)r2 (20) 


In this method of initialization, you set P and P’ near r=0 and then integrate 
step-wise out to large r. At large r, if the energy, E, you chose was an allowed 
energy (an eigenvalue), the reduced wave function, P, will go smoothly to 
zero as will, therefore, the radial wave function R(r) = P(r)/r. (This is the 
behavior of P and P’ whenever V(r) is dominated by the angular momentum 
term 2(L+1)/(2r2) for small r. The asymptotic equation for small r is then 


Pp’ = aA P which has the solution P = pets) 


A SCREENED COULOMB POTENTIAL PROGRAM 


* One implementation of this computer strategy, based on initialization by 
expansion for small r, is shown below. The potential is a screened Coulomb 
potential which is relatively weak; the screening distance is 8 Bohr radii = 4 A. 


GET-3DSCHe 
LIST 
3DSCHe 


10 DEF FNVCR) =-ExP(-R/a>/R—- Potential, Vir) 
20 PRINT “RCSTART),RCEND)» DELTA*R, ORBITAL @#(L)>,» RADIAL #(N), ENERGY?" 
38 INPUT RBsR9,R7,LENGE 

42 LET PO®EXPC-RB/N)D REI CLs £) 

SO LET Pl=POe(-RB/NO(L+1)/RB) > (nitialization at small r 
6@ LET RRB 

78 PRINT Rs PGs PS/P 

8@ LET RS@(R9-RB) /26 

98 LET R6@ReR7/2 


198 LET Pes2e( FNVC RO) +L *(L4 1) /( 2 R6eR6) -E)*CPO+P1eR7I/2) | Integration of 


11@ LET PGePG+(P1+P2eR7/2)R7 the Schrodinger 
128 LET PisPieP20R7 ( ; 
136 LET R®R+R7 ’ equation 


14@ IF R«RS THEN 170 

15@ LET RS@R5¢(R9-R8)/20 ; Print loop 
168 PRINT RPO, PA/R 

170 IF R<ROI THEN 96-Aeturn for next Ar 
182 PRINT 

198 GOTO 28 

20@ END 
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The variables PO, P1, and P2 are the values of P, P’, and P’ for each radial dis- 
tance r in turn. The print-test is not a necessary part of the strategy, but it saves 
printing out every answer—all several hundred of them for each energy guess. 


The ground state is expected to have the orbital angular momentum quantum 
number, & = 0. If some other state has lower energy, you'll find out when you 
search for that state. The unscreened Coulomb potential (the hydrogen atom) 
has energies (in atomic units) E = -1/(2n2). So we hunt near E = -.5; the dif- 
ference from E = -.5 that you will see in a moment is real and is due to the 
screening. 


You need not use the analytic result for the hydrogen atom to guess your 
energies; you could look at experimental results instead. From optical experi- 
ments on the spectrum of hydrogen, the energy levels can be written as 
E = -Ryd/n2, where Ryd = .5. You need not have a good guess for your first 
energy trial; you'll home in on the eigenvalue sooner or later. It helps to think 
about your first energy guess because you will home in on the answer 
more quickly. 


SAMPLE EXERCISE 
Let’s walk through a typical exercise. For the potential V(r) =-exp(-r/8)/r: 


1. Plot the 2=0 (s state) and 2=7 (p state) effective potentials from r=0 
through r= 5, 


The program is very simple. A printing version might look like this: 


3DPOT 


16 FOR X®@=#.-2 TO 5 STEP -2 
28 LET Ve=-EXP(@X8/8) /X8 

38 PRINT X@> Vs V42/( 2eXB#XB) 
48 NEXT X@ 


58 END 
END 
RUN 
3DPOT 
e2 ~4e87655 28.1235 
6 -~12-54624 1223154 
08 -1.2-13165 431453 
I 2882497 117583 
le2 -2717257 -2e28 122E-82 
164 “0599612 “8694808 2E-82 
1-6 -.5117087 -~-12!1¢0@82 
1-8 ~ + 44362 "« 134978 
2 2 3894 -¢ 1394 
202 "se 34526 =e 138648 
2e4 ~- 308674 -2135@63 
226 “277835 -- 1299 66 
228 ~2251674 -. 124123 
3e ~2 229897 “2117985 
3e2 ~~. 289475 -2111819 
324 ~. 192285 -~- 18578 
36 ~e I T7119 -9 99S85E-B2 
3-8 “+ 163654 ~9.- @4BIDE-B2 
ae “0 151633 -8-91327E~B2 
4e2 ~« 1408 a7 -6-« 41IS5S72E~82 
aed “e 131125 -T-9ATZIE- 82 
426 ~e 122327 -7.S58682E-@2 
4.8 #2114336 -» 870933 
Se “2187852 ~6278523E-B2 
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Terminal plotting (as discussed in the Appendix) could also be used. The 
changes in the program and output might look like: 


S$ READ X8sX9eY8LY9 

6 DATA 8.527141 

7 GOSUB 99886 

18 FOR X9e.2 TO 5 STEP .2 
26 LET Ve-EXP(-X68/8) /X8B 
36 LET Z8#1 

31 LET Y@=2V 

32 GOSUB 9198 

33 LET Zg#2 

34 LET YOu Ve2/( 2eXGeXO) 
35 GOSUB 9188 

48 NEXT X@ 

4l LET Z268Xde7 G6 

42 GOSUB 91868 

45 GOSUB 92890 

$8 END 


APPEN D- TTY PLO 
RUN 
SDPOT 


MAX Y= i 
RXAXKXXAXAAAAARAANKA KX AXKXXXAKAAKAAKAHAA AAR KKKM MK KAHAN AK KA AK NK NHK KK KKK KKK 


CRAMMER KKK EH 
to 

© 

© 
~ 1 
~- 1 
—~ © 
- © 
- 9 
-— 


XXXXXXAXAXKXAKAAKARKAKAAK AXA A KK XXX XXX XX XAKKK KAA XKKXAKAKK KKK HHH KNX KX KKK KKH HK 
MIN Y¥=2z-} MIN Xe 6 . MAX Xe 5 


We use the second three-dimensional Schrodinger equation program dis- 
cussed above (the one which initializes P and P’ near r=0. The energy of 
the hydrogen atom for this state is E = -.5. From the RUN you can see 
that the equivalent screened Coulomb energy is about -.384. So the 
state has been shifted in energy due to the screening. The wave function 
is very similar to that of the hydrogen atom at small r but falls off faster 
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2{(a} Sketch the radial wave function of the lowest =O state. Now find the 
energy and wave function for this state with the computer. 


than the unscreened wave function at large r. 


RUN 
3DSCH2 


RC START). RC END)» DEL TA-R, ORBI TAL 


TLE-8s 165 681.80 1lse-05 


}- GBBABE=88 
5] 

1-281 
1-51 
2091 
2-59981 
3+ 88882 
3-58803 
4e68004 
42 58686 
Se 98827 
5-52808 
6-6 38889 
6.5201 
72¢8881) 
7258812 
Be-8GGia 
8-S5@B15 
9.-88216 
9258617 
$@28002 


RC START)» RC END)» DEL TAR, ORB! TAL 


1+ @BBBGE~88 
°31G115 
o 398447 
«392198 
«38 1859 
391614 
244188 1 
2552675 
752822 
1» 88464 
1262346 
20487248 
386882 
6.808183 
926315 
1$- 3447 
24-5615 
392 4649 
63¢ 6137 
182-815 
166.556 


PLE*8, 10228ls Bs lame 38 


1-08688E-28 
05] 
1-81 
165! 
2e21 
2258881 
3- 28882 
3e 58963 
4e B®8884 
4e S8OBE 
Se 08887 
Se 58008 
be 22889 
60568! 
7e@OG811 
7+58012 
8- 28814 
8+58815 
9228816 
936817 
14-6 @882 


RCSTART). RCEND)» DEL TA-R, ORBI TAL 


1-8@8QGE-88 
3985322 
2 3624877 
319217 
© 241372 
« 159243 
TeP2664E-Z2 
Je 6999 7TE“83 
721799 3E-@2 
me 154643 
-« 255851 
-+2 38 6882 
-2569 487 
2838197 
#1626895 
-1-76683 
-2.59 68 1 
~3-8 1461 
-5- 63528 
~8.935335 
~122e4177 


VIE-8. 18,80-01,02 14°0384 


1. 88S8GE-28 
31 

le @! 
1.51 
2-81 
2258821 
3e #8802 
3 58883 
4-28064 
4e 508866 
$e 88007 
5-58608 
6e 88089 
6- 5881 
7¢@881t1 
7¢58812 
8-80014 
856815 
9-88616 
9.58817 
16-9862 


R¢ START)» RC END) » DEL TA~R, ORBI TAL 


? 
END 


le @GBBGE-25 
© 386485 
+ 36873 
0335219 
271278 
«28678 
2149724 
© 183993 
6e88$428E-82 
4. {B3S9E-2 
1-88 1L1ISE-82 
“8. IBI61E-84 
-2-B85754E-82 
~4.95395E-¢2 
-74 348 S4E-22 
~2e 113534 
-2 177838 
-+ 268917 
-. 487663 
-2 618411 
0939 641 


le 

« 688869 
- 386582 
2° 259734 
2 18998 
° 156645 
2 1a7293 
« 157986 
« 186804 
«241829 
« 324689 
© 644793 
935528 
1637591 
2.804593 
3e 87814 
4.64285 
7 - 86887 
ta@.8225 
16-6553 


le 
° 59867 
- 358888 
211482 
«1288686 
6-369 68E-22 
» 826422 
le @5713E-23 
-1. 79949 6E-Be2 
~3+43686E-82 
-5. 182809 6E-682 
-~7e@3412E-2e2 
+9 - A899 7E- B2 
~e l27721 
-.172784 
0235574 
-2 3238 46 
=« 4487 
-2 626131 
-«879 285 
-1-24175 


le 

» 686794 

¢ 365879 
2221999 

« 1349 64 
Be271L1BE-B2 
4099 B77E-22 
229712 

lo FISSSE-B2 
Ge LIBSBE-83 
Je 76239E~- 83 
“Le &73I36E-84 
-3- 429 1BE-83 
~60 698 27TE*@3 
~12¢84977E~O2 
-125a842E-82 
~20¢21294E- 62 
#3. 16367E-62 
~4-52951E-82 
= 6¢ 589 ABE-82 
“9239 624E-82 


#(L.). RADIAL 


#¢L), RADI AL 


@¢L).s RADIAL 


6(L), RADI AL 


@(N), ENERGY? 


#(N), ENERGY? 


#(N)s ENERGY? 


@(N)s ENERGY? 
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2.(b) Sketch the radial wave function for the first excited °=O state. Now find 
the energy and wave function for this state with the computer. 


We use the same program. The energy for a hydrogen atom in this state 
is E = ~.125. From the RUN you can see that this state is almost lost to 
the continuum (i.e., E = -.0375 ts nearly positive). 


RUM 
IDSCH2 


R¢ START)» RCEND),» DELTA-R, ORBITAL #(L)» RADI AL 
2 1EW8 4 18s 081285 22-0 125 


1-8889GK-88 
«$1 
1-81 
1-5! 
2081 

2- 569981 
3- 88602 
Je 56883 
4- 68964 
4.58086 
§- 88067 
5.599008 
6- 88889 
6.5601 
7-8G011 
7+S$8812 
8.86014 
850815 
Ge BBG16 
9250817 
10-8802 


1-S6800E-08 
° 298 287 
+ 322958 
+ 228945 
60979 ISE~G2 
~8-9OG2S5E-82 
=. 244526 
-° 387851 
0 $19923 
-» 645218 
-2776876 
~e9 8191 
- 1284696 
-1.22838 
~ 12426356 
-1-¢67963 
~ 1299489 
-2+387 68 
~2e88 243 
-3.58595 
~ 4. 29388 


le 
38472 
319167 
e 14632! 
3- AT232E~82 
~3« 56889 E-92 
~. 88 15908 
-2e118614 
~2 129979 
=« 14338 
~e 1§4813 
~e 163981 
~e 174824 
~. 187748 
-2 283792 
~e 2239 47 
~e 249287 
-< 288899 
~« 328265 
--3000a 
=.+ 4293 


RC START). RC END)» DELTA-R, ORBITAL #(L),. RADI AL 
TLE-8, 200-012 0s 22-283 


1- 8SGSGE-88 
1-8) 
2-81 
3-89802 
Ae BOBO4 
5288887 
62 88809 
7¢808i1 
8-O8GG14 
9-B8GG16 
18.8692 
11-8002 
1228882 
13-0882 
14-8883 
15-8963 
16e 88083 
17-690! 
18.81 
19.-0@98 
28-0896 


RC START) » RC END) » DEL TA~Rs ORBITAL 


le OBABBE- BE 
¢ 386865 
1+ 8S067E-62 
~. 328202 
-° $3714 
-+ 626289 
-- 6190881 
-.552657 
~» 44896 
=» 332253 
“2211838 
“928869 1E~82 
2e98BS37E-82 
» {3861 
© 2759589 
#411125 
° 562671 
«738 602 
*9569 62 
le 20758 
12652516 


2 1LE"8s 26s 6 Bis Bs 22-2 B38 


le BBEEBE- BE 
Le 81 
2081 

Je 88682 
4.89894 
Se 88807 
6- 88889 
7-860! 
8.98814 
9-88816 
18-8882 
11-8862 
12-8882 
13-@8¢2 
14.8883 
15-9883 
16-8883 
17-8981 
18-@1 
19-8898 
20-3896 


RC START) » RC END)» DELTA~R» ORBI TAL 


? 
END 


1e BBGS9BE-B8 
©3969 13 
1-29796E-82 
-+ 316965 
-« 537928 
~« 636199 
~¢ 64284 
=-.592335 
-e512354 
=» 428293 
-+ 32689 ] 
-«237352 
-° 153334 
~7+42386E-82 
2eO8G354E-83 
7+839B4E-82 
«158723 
e 247585 
«351452 
© 47562 
« 627367 


be 
#3836395 
40978 45E-83 
*2° 186733 
72134264 
“2 l25256 
-~e 183165 
-+ 878664 
-.856119 
“3+ 69 1 63E-82 
~2- 1I@34E-G2 
“8. IS7T9SE-“93 
2+ 48776E-83 
le 1S852E-82 
219711 
2+ 74878 E-82 
3¢51663E-82 
4. 34468 E- 82 
9* 28 819E-82 
6-6 352a2E-82 
7+ 62213E-82 


le 
¢ 3838 75 
6e A57S3E-83 
-<185654 
~0 13448 
=e £27238 
me 187138 
~8- 46466E-82 
= 6¢ 48432E-82 
~Ae 6698 AE~ 82 
-3+ 2688 GE-82 
~2e1S77TILE-82 
~1e27776E-62 
“5. 7899 4E-83 
be A31G7TE-84 
S« 2259 2E-83 
9-928B2E-83 
le 45619 E-62 
1-9S143E-82 
2- 4988 1E-62 
3e 13532E-82 


#(L). RADI AL 


#(N)s ENERGY? 


#@(N). ENERGY? 


@(L)sRADIAL #{N),. ENERGY? 


#(€N), ENERGY? 
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Plots of these two 2=0 states can be made on the terminal (as discussed 
in the Appendix). The changes in your program and a RUN are: 


§ READ X8sX9sY8s YP 
6 DATA 62 185-42. 1 
7 GOSUB 9080 

70 LET X9#R 

71 LET Y9=P6/R 

72 LET 2767! 

73 GOSUB 9182 

166 LET X@=R 

161 LET Y@=P0/R 
162 GOSUB 918d 

175 LET Z@*X@aY9«8 
176 GOSUB 9186 

188 GOSUB 9268 
APPEN D- TTY PLO 

RUN 

I3IDSCK2 


RC START)» RCEND) » DELTA=Re ORBITAL #(L).s RADI AL #@(N)», ENERGY? 
PIE“6, 105001. Os lave 384 
MAX Y=? 1 
XXXXMKXMKXKK XKXKX KX KKK K XKKK KKM KXKKK KX KKK K ARKH XKHKK MAK KK KKK HK KKKK KX XK XKXKK X 


ROS KKK RK KKK KKK KKK KKK KKK RK KK 


RXXXKXKXKXAXKKKAAKRAAKRAK KA ARK XK KK AKK KKAAKKAKMK KAR KK AAK KK KK KKAKK KKK KR AKAKK K AKAKK AKKK 
MIN Yer.2 MIN X= 6 MAX X= 18 


RC START) » RC END)» DELTA-Rs ORBITAL @(L)s RADIAL #( i). EXERGY? 
218-6. 162 681s Os 22-0835 


MAX Ye 1 
AK AKA AAAs sos sro teste se sesrsososoedrtrtrsrctide trode drtath Arras sadodra Pate Potodtodied * o.¢ Oot 


hq Hh Hd Hh He eh He Hh ee Ke eh oe ee ee eee ee Ee 


AXXRAXXAXXKXAXXAKAKAXKAAKAK AKA ARK AKA AKAKAKAAKAAK RA RAK KAKA AK KKK AANAAK RAK AX AXAAKA HK 
MIN Y=s-.2 MIN X= @ MAX Xe 196 

RC START)» RC END)» DELTA“R, ORBITAL @(L)>s RADIAL #@(N).s ENERGY? 

? 

END 
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2(c) Sketch the radial wave function for the lowest energy 2=1 state. Now 
find the energy wave function for this state with the computer. 


We use the same program. The hydrogen atom energy for this state is 

= -,125. The lowest p state and the second lowest s state are degenerate 
in hydrogen. This state is again nearly lost to the continuum for this 
screened Coulomb case. If the state were lost to the continuum, a RUN 
for E=0 would not produce enough curvature (the kinetic energy term 
in the Schrodinger equation) to bring the tail to zero. This state starts 
out from R(0)=0 and has one (positive) lobe before the tail goes 
smoothly to zero. The RUN shows that E is about ~.0325 atomic units. 


RUN 
3DSCHe 


RC START) » RC END)» DEL TA~R, ORBI TAL 


TLE-8,26.+81s122690 


1. @8886E-28 le SOGBGE- 16 1+ BBSBBE- BE 
te 81 Je 847 30E-86 Je BIL7FISE-86 
2-81 7226512E-96 de 61449 E- 86 
3e 83882 9+ 636 28E-86 3eLI2TAE-86 
&e 80804 9 e83I2ZAE-96 2e 488 26 E-96 
Se 88007 $e 26 37GE-66 1e 6567 2E-86 
62 89889 Se G1 AG2E-86 9-2 33565S5E-87 
7069611 2-389 A2E-86 3e 3788 4E-87 
$e 888i4 “1LeBIIIE-66 163641 1E-07 
9289816 ~ Ac 48 AA6E-06 ~4c98265E-87 
18-8002 -7- 6769 BE~@6 =7- 67876E°-97 
11¢8882 =1¢85752E-95 -9+61363E-°87 
12-6662 -1¢31766E-65 -1¢89 7S3E-86 
13-8862 ~1-S460S5E-95 =1o1898SE-86 
14. 8983 -)]-746S8E-85 = $e 2Aa7S3E-86 
15-8863 -1-9216GE-85 -1- 28 164E-86 
16-6683 2087441 E-95 -}e 29 648 E-86 
17-8961 ~2e 268 32E-85 ~1- 299 88E-86 
16<8! ~2¢ 32757E-95 ~1e 89298 E-86 
19-8098 ~2< 43268 E-65 -~1-279 7T8E-86 
26-809 6 -2¢52745E-85 ~1-26312E-86 


#(L).s RADIAL @¢(N)s ENERGY? 


RC START) » RC END)» DELTA“Re ORBITAL #(L). RADIAL #(N)s ENERGY? 
? 
END 


Since these last two states have different energies, the degeneracy of 
the unscreened Coulomb potential is broken by even this fairly weak 
screening. 


3. Calculate the difference in energy between the two X=0 states in atomic 
units, eV, Hertz, and A. 


(1 eV = 1.602 x 10-12 ergs; h = 6.624 x 10°27 erg-sec; c = 2.998 x 1019 
cm/sec) 


E = - .038 + .384 = .346 atomic units = 9.41 eV = 1.51 x 10°"? ergs. 
frequency = 2.28 x 1015 Hertz; lambda = 1317 A. 


Note that the wavelength for emission from this transition compares 
well with the Lyman alpha line (1216 A). The difference between the 
second &=0 state and the lowest £=1 state is .005 atomic units=.136 eV. 
The spectral line for this transition has lambda=91110 A (in the infrared). 


4. Plot the probability distribution of the full wave function w{r,é,o) = 
R(r)¥ (0,0) in the 6=0- $= 180 plane for the ground state (the lowest-lying 
R=0 state). 


Terminal plotting can be used. The following program stores the radial 
and angular parts of the wave function and then finds the current com- 
bination for each point plotted. 
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3DR-T 


19 

28 

30 

48 

52 

68 

70 

88 

98 

188 
118 
128 
136 
148 
158 
168 
178 
186 
198 
286 
218 
228 
230 
248 
258 
268 
276 
280 
298 
388 
318 
328 
338 
348 
338 
368 
376 
388 
396 
488 
418 
426 
436 
446 
456 
466 
A7@ 
488 
492 
598 
516 
3268 
338 
548 


DEF FNVC(R) ©-EXP(-R/8) /R— Potential, V(r) 

DIM TL401}.RC 402) 

PRINT “RCSTART).»RCEND), ORBITAL #(L)>. ENERGY ?” 
INPUT RB, R9oLZE ) 

LET R7=( RO-RS8) /ae seb atd ee. oie 

Let t7=3.14189/480 ( /Nitialization 

LET RO=To=6 \ 

FOR I=! TO 491 

LET Ta-30 14159 /2+¢1-1)"T7 

LET TC1I=1 Calculation of 
IF L=@ THEN 136 Q)- 

LET TC1iJ#( COSC T))*2 ©(8), : 

IF TC1)<T@ THEN 158 stored in T{ ) 
LET T@#TC1) 

NEXT I 

LET R@RE 

LET P@#Rt(L+1) a ae 

LET P1=(L¢1)#RIL Initialization of R(r) 
LET I=1 

LET RC 1)" POeP8/(R*R) . 
LET R6@R+R7/2 Integration 
LET P2=2e(FNV(R6)+L*(L+ 1) /( 2#R6eR6)-E) *¢Pa+Pi*«R7/2) | Of the 

LET P@@PO+ (P14 P2eR7/2)*R7 


LET PisP1+P2#R7 | Schrodinger 
LET R=R+R7 ( R2(r) equation 
LET late! 


LET RC1}=Pa*Po/¢ReR)! 
IF RCII<RO THEN 380 | 

LET RO@REI) - Find [R2(r)/r2] 
IF R<«RO+R7/2 THEN 216 | 

PRINT "MAX. Re"3ROS" MAX. Tas TS 
PRINT “LEFT X» RIGHT X» BOTTOM Y. TOP Y. MAX. FOR 20?" ‘eiafj 
INPUT X8.X9sY8s 9229 | Initialize 
GOSUB 9200 | plot 

FOR X@*X8 TO X9 STEP (x9-x8)/14— Step across x 

FOR Y@=Y8 TO Y9 STEP (Y9-Y8)/18— Step across y 

IF ABS(X@)+ABS(Y@)=@ THEN 490 

LET R=SQR(XG*XG+Y B*Y BD 

1F R>CRO-R8) THEN 49@ 


Store ——z—in R{ } 
r 


max 


LET [@@1¢INTCR/RI+.S) Find right elements 
IF X® <> @ THEN 449 

LET 1101 of R( )andT{ } 
GOTO 450 


LET 11m2Q1+INTCATNCYO/X@) /T7465) 

LET ZOsINT(9#RC103#TC11)/29¢.5)— Find right character for |p|? 
IF 288(Z@-8) <= @ THEN 48¢ 2 

ee cece } Force all too-big |W |? to be 9 
GosuB 9186-—Store point in plot 

NEXT Y@ 

NEXT X@ 

gosus 92ee—~Print out the plot 

PRINT 

GoTo 328— Return for new plot 

STOP 
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RUN 
3DR-T 


RC START)» RCEND)»sORBITAL #C€L).» ENERGY 7 

7P1E-8. Ss Bo -+« 384 

MAXe R=® 9758235 MAXe T# 1 

LEFT X, RIGHT X,. BOTTOM Ys TOP Ys MAX- FOR 26? 

Pele As be Ass lesl 

MAX Y= 1 
XXXXXMAKXKAKKKKKKXKAKAKK RK KKAK KK KK KKM KARR RK KK KAKK KKK KKK KX RKHKM KKK HK KKK KK RIKK 

Yo 8 i 1 1 I I I i J j l 1 @ 
Y 

Y 

¥@ l 1 \ I i 2 2 2 j i 1 ! 1 
Y 

Y 

6 1 1 I 2 2 3 3 3 2 2 1 1 1 
Y 

Y 

Y@ i i 1 2 3 4 4 4 3 2 1 l 1 
Y 

Y 

Yt 1 i 2 3 4 5 6 5 4 3 2 I 1 
Y 

Y 

Y! 1 I 2 3 4 6 9 6 4 3 2 1 1 
Y 

Y 

Yi 1 l 2 3 4 5 6 5 4 3 2 1 1 
Y 

Y 

YQ 1 1 l 2 3 4 4 4 3 2 1 I 1 
Y 

Y 

YS i } 1 2 2 3 3 3 2 2 1 } 1 
Y 

v 

Y@ i I 1 1 1 2 2 2 i 1 1 I ! 
Y 

Y 

Yé a 1 1 1 1 J I 1 1 1 ! 
RXXXAXXXAAKAAXAXAXAAKAXAXAXKAAAAAKAKAKAAKAAAA KAKA MAMMA KKK HN KM AMM MAM K HMA MK RAH NH HM 
MIN Ye-1i MIN X#-1].4 MAX X*® 1-4 

LEFT X» RIGHT X». BOTTOM Y,. TOP Ys MAXe FOR ZG? 

? 


END 


EXERCISES — Use computer when appropriate. 
|. For the screened Coulomb potential V(r) = -exp(-r/10)/r: 
1. Plot the effective potential for 2= 0 and 2 = 7. 


2.{a) Find the radial wave function and the energy for the lowest energy 
2 = O state. 


2.(b) Find the radial wave function and the energy for the first excited 
2 = Ostate 


2.{c) Find the radial wave function and the energy for the lowest energy 
2 = 7 state. 


3. Calculate the difference in energy between the two %=0 states in 
atomic units, eV, Hertz, and A. 


4. Plot the probability distribution of the full wave function (r,6,¢) = 
Ri(r)Y¥ (0,0) in the ¢ = 0 - ¢ = 180° plane for the lowest lying p state. 


Hint: Use the program which initializes P(r) near the origin. 
ll. For the screened Coulomb potential V(r) = -exp(-r/5)/r: 
1. Plot the effective potential for = O and = 7. 


2.(a) Find the radial wave function and the energy of the lowest 2=0 
State. 
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2.(b) Find the radial wave function and the energy for the first excited 
L = O state. 


2.(c) Find the radial wave function and the energy for the lowest energy 
R= 7 state. 


3. Calculate the difference in energy between the two 2=0 states in 
atomic units, eV, Hertz, and A. 


4, Plot the probability distribution of the full wave function W(r,6,) = 
Rir)/Y(6,0) in the ¢=O0-o= 180° plane for the lowest-lying p state. 


Hint: Use the program which initializes P(r) near the origin. 


_ For the Lennard-Jones 6-12 potential, V(r) =400 (1/r'2- 1/r®): 


1. Plot the effective potential for 2=O0 and &=7. 


2.(a) Find the radial wave function and the energy of the lowest 2=0 
State. 


2.(b) Find the radial wave function and the energy for the first excited 
R= 0 state. 


2.(c) Find the radial wave function and the energy for the lowest energy 
R= 7 state. 


3. Calculate the difference in energy between the two =O states in 
atomic units, eV, Hertz, and A. 


4. Plot the probability distribution of the full wave function W(r,0,o) = 
Rir)Y¥ (0,6) in the 6=0- o= 180° plane for the lowest-lying p state. 


Hint: Use the program which initializes P(r) at large r. 


iV. For the screened Coulomb potential V(r) = -exp(-r)/r: 


1. Plot the effective potential for 2=O and = 7. 


2.(a) Find the radial wave function and the energy of the lowest 2=0 
State. 


2. (b) Find the radial wave function and the energy for the first excited 
L=O state. 


2.(c) Find the radial wave function and the energy for the lowest energy 
R= 7 state. 


3. Calculate the difference in energy between the two =O states in 
atomic units, eV, Hertz, and A. 


4, Plot the probability distribution of the full wave function w(r,0,6) = 
Rir)Y (8,0) in the 6=0- = 180° plane for the lowest-lying p state. 


Hint: Use the program which initializes P(r) near the origin. 


*V. For the potential family, V(r) = 400 * (r-12-r-™): 


1. Find theground state energy and wave function for m=3, 4, 6, and 10. 
2. Plot the energy of the ground state as a function of m. 


3. Find the small r and large r behaviors of each wave function. 
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4. Can you extrapolate your results to all m? 


. For the potential family, V(r) = -exp(-r/ry)/r: 


1. Find the ground state energies and wave functions for r=1, 2, 4, 8, 
and 12, 


2. Plot the energies versus the screening parameter, ro: 
3. Find the small r and large r behavior of each wave function. 


4, Can you generalize your results to all screening parameters, lo? 


**VII. For the potential family, V(r) =-1/r™: 


1. Find the ground state energies and wave functions for m=1, 3, 6, 
and 10. 


2. Plot the energies of the ground states as a function of m. 
3. Find the small r and large r behavior of each wave function. 


4. Can you extrapolate your results to the behavior of these ground 
states for all m? 
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CHAPTER THREE: THE ONE-DIMENSIONAL, TIME-DEPENDNT, 
SCHRODINGER EQUATION FOR TIME-INDEPENDENT POTENTIALS 


INTRODUCTION 


The time-dependent Schrodinger equation in one dimension ts 


h2 a2y | my 
—-—— — + = — 
am ax2 Vi(x,t) w if at (21) 
If the potential, V(x,t), is independent of time, V=V(x), then the equation 
separates by setting (x,t) = P(x)T(t). The equation for P(x) is the time- 
independent Schrodinger equation 


h2 d2P 
The equation for T(t) is 
dT <E 
-rohaiad ~ T (23) 


E is the separation constant and, as you can see from its position in Equation 
(22), plays the role of the energy. 


The equation for T(t) (Equation (23)) solves immediately to 


=t 


T(t) = Ae! (24) 


which in Hartree atomic units (see the Preface) is 


T= Ae" (25) 
The constant A is just a normalization constant. 
The spatial equation (Equation (1)) can be solved analytically (for a few 
problems) or numerically as discussed in Chapter One. This chapter empha- 
sizes the effects of the time part of the full wave function W(x,t) = P(x)exp 
(-iEt). 
STATES OF THE INFINITE SQUARE WELL 
The infinite square well, 
co «6X <0 
V(x) =" O O<x<L (26) 
co XSL 


can be solved analytically. The result is standing waves whose spatial wave 
functions are 


P(x) = B, sin{nax/L) n= 1, 2,--- (27) 


B,, being a normalization constant. The full wave function for these discrete 
states are then 


v,, (x,t) = C.sin(nax/L) e'Ent (28) 
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The energies of the states are 


_ 1f7¥ 2 
E, = 3(2) 0 


THE BEHAVIOR OF SETS OF DISCRETE STATES 


n=1,2,--- (29) 


Consider setting upa particular shape for the electron’s distribution at time = 0. 

For definiteness, let’s set up (an approximation to) a square wave with 

% wavelength in the well. By the ideas of Fourier Series, we know that we 

want to add up a number of these discrete states (in particular, those with 
= 1, 3,5, 7, ...) with varying amounts for each n. If we write 


square wave = 2A, P, (x) 


then theA, ‘s are proportional to 1/n in order to produce a square wave. 


After setting up such a shape at t=O, you then let time start passing. Each 
mode of the well will oscillate in time with its own angular frequency w, = 
E./h or w, = E, in atomic units. The energy depends on the wavelength 
squared or the wave vector of the state squared. The higher energy modes run 
ahead of the lower energy modes, and the shape distorts away from a square 
wave. 


This result is unlike the result of a square standing wave on a stretched string. 
On a stretched string, the speeds of the various waves are independent of the 
wavelengths (or wave vectors k, = n7/L for the discrete states). The speed, v, 
of waves on a string depends only on the parameters of the string (tension, etc.). 


In the case of these electron waves, the speed of the wave, v = w/k = E/hk = 
hk/2m or v = k/2 in atomic units. The speed depends on the wave vector 
(that is, the wavelength). This is an example of what is called a dispersive 
medium. The vacuum (where V(x) = O everywhere) is a dispersive medium 
for electron waves. 


Let’s watch the time development of our electron square wave in the infinite 
square well. A program (and its RUN) to perform the sums and products 
necessary might be: 


WELSET 

t@ DIM AC1S,21),EC15),CC15),S015) 
26 LET P1=3-14159 

3@ LET Lei Length of well 


a@ FOR I#1 TO 18<—Step through 15 states 
$6 LET Neoe(l-1)41 

68 LET KO@eN#PI1/AL . 

76 LET EC13=Ke«Ke/2«ENergy of each state 


8@ FOR J#1 TO el Step across X 
96 
t 18@ 


Store energies 
and spatial 


ET XO=(J-1)*L/28 é 
LET X@=(J wave function 


LET ACI, JISSINCKO8*XB) /N-«*- Snatial wave function 


116 NEXT J 





122 
13@ 
140 


156 
168 
( 178 
188 


—* 198 
288 
2108 


228 
238 
2468 
258 


268 
276 
288 
298 
388 


NEXT I 

PRINT “TIME?” 

INPUT TO 

FOR I#i TO 15 

LET ClLIJ@COSCECII*TO) ‘ 
LET DSc cel i sta Store exp (-iEt) 
NEXT I 

FOR Jel TO 21<«—Step across X 

LET X@e(J~1)*L/2@ 

LET AL@A2@@ _ Step through 15 states 
FOR l#! TO 15 

LET AlwAl+ALI,U)*Cr1)-<«— Real part of W(x,t) 
LET AZ*A2+ALI2J)4*SC1 , 

sh A2+AC ~~ Imaginary part 
PRINT XQ, ALPAL+A2KA2 

NEXT J 

PRINT 


ais 138 Return for new time 


For each state 


RUN 
WEL SET 


TIME? 


085 
a} 
015 
2 
25 
3 
35 


2 45 
25 
«55 
«6 
2 65 


75 


25 


35 
2& 
245 
25 
55 
26 
«65 
07 
75 
8 
85 
«9 
095 


TIME? 
72 


85 
e} 
015 
22 
25 
3 
35 
«4 
«45 
oe) 
288 
06 
« 65 
0? 
75 
«8 
85 
09 
«95 


TIME? 


END 


8 
«646566 
782685 
«613192 
¢ $7 3336 
618871 
« 649561 
» 616355 
e $69 667 
© 616991 
© 643279 
«616989 
© $89 667 
e 616355 
« 649561 
618872 
$7 3335 
61319 
782685 
« 646569 
16 86422E-12 


8 

«189818 
6-219 48E-82 
© 295673 
$5187 1 

e 698341 

2 434285 

° 689529 
©7489 19 
1-38 383 
1-88855 
1.58383 
74892 

« 689528 

© 434284 
«69634 
«55187 

« 298672 

6- 21952E-92 
«199818 

1-7 4668E-13 


8 

4.38 285E-82 
9-8189 6E-82 
« 3445 67 

+ 148464 
651564 
942334 
1-23172 
}-3574] 

2+- 25539 
$2357.41 
1e23172 
942334 
65156! 
«148464 

+» 344566 

9-8 1895E-82 
4-38 28 6E-82 
e 124333 

5- 16428E-13 
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Terminal plotting (as discussed in the Appendix) can also be used. The changes 


you might make in the program and a RUN follow: 


5 READ X8.X92Y8sY9 
6 DATA 0.1,0,2-5 

7 GOSUB 9806 

268 LET YB#AIFAI+A2#AC | 
278 LET 2821 

276 NEXT J 

262 LET Z@#1 

264 GOSUB 91886 

288 GOSUB 9268 
APPEND-TTY PLO 

RUN 

WELSET 


TIME? 

28 

MAX Y® 2-5 

XXXXXKXK AKAN XXX XAXAK AKA KAKA ARK KX AX KX KX XAXXAKAAK KA AKA KAK KKK AKAN K KKK AK AKK K RKX A AKX 


MRK RK KKK KEK KKK KKK 


1 i 
XXXXXXRXXXAAK ARAN KAKA HAKKAR AAARAAAKAAARAAAMAK AAA KAMA MAA A MMA MAM MMMM 
MIN Ye" 6 MIN X@# 6 MAX X# | 


MAX Y= 2-5 
n4,8.5.9.9,0.9.9.0:9 9.0.9.9.0.0. 09.9.0. 90.0 9.9.0 0.9.6.60.000.0.0.0:9.0.09.9.0:9:9.0.9.90.9:9.19 9.0.0 9009,99,9,9,0, 0,09 0.9.9 


ae | } i 
i i 
XAXAAAKRKAANRAAA NAKA ARAMA AAK MA KAAARAA KH ARAAAAKRAKAAAAKAAK AMAA MM MAA MAMMA MAA 
MIN Y= @ MIN Xe 6 MAX X= I 


RRR EK KKK KKK EK KEK KK KKK KKK KK 
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TIME? 


Max Y= 2.5 
XXXKAXAAXARAKARAAAAKA KAKA KAKA AAAAAAAKAAKRAMARARAAAKAAAAANKAAAAANM MAMMA GAMA 


SEK eK KKK KEK KK KKK KKK KKK 


t 2 1 io! ] 
YI i 
XAXXXXXXAXKK AKA XAAX AAR AK RAK AKAN XARA KAKA AKAK RK X KKK KXAKA KK KK KKK KAMAKAK KAN XX AXA XAARK 
MIN Ye @ MIN X= @ MAX Xs 1 
TIME? 
? 


END 


This same sort of procedure can be carried out for the discrete states of any 
one-dimensional potential. If an analytic solution to the problem exists, it Is 
usually more efficient to use that. If an analytic solution does not exist, you 
can solve numerically for the spatial wave function, P(x), as discussed in 
Chapter One. The same sort of procedure can be carried out for states in a 
continuum (for example, the states of a finite square well whose energies lie 
above the top of the well). Before we discuss the behavior of sets of continuum 
states as time passes, we need to discuss the spatial wave function solutions 
for one-dimensional continuum states. 


ONE-DIMENSIONAL CONTINUUM STATES 


Suppose we have some potential well such that, for |x| > xp), the potential, 
V(x), is essentially zero. Then, for | x | > xg, the Schrodinger equation for the 
spatial wave function, P(x), becomes 


P(x) = -2EP(x) (30) 


The solution for P(x) is P(x) = exp(+ikx), exp(-ikx) or, equivalently, P(x) = 
cos(kx), sin(kx). Neither solution must be discarded because of boundary 
conditions: both solutions are well-behaved everywhere; there is no classi- 
cally forbidden region for large x which forces one or the other solution to be 
unphysical. There are two solutions to the equation for every energy, E. 
Furthermore, since there are no boundary conditions which limit the accept- 
able energies, every energy, E, is allowed. Hence the states form a continuum 
of pairs of states. 


\f the potential, V(x), is symmetric, the theorem discussed in Chapter One 
(and the Appendix) tells us we can choose the states to be wholly even in x or 
wholly odd in x if we wish. Since any linear combination of solutions to the 
Schrodinger equation is itself a solution, we could choose any arbitrary phase 
and make sin(kx+@) one of our two solutions; the other solution (by ortho- 
gonality) would then be cos(kx+@). 
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The simplest potential of the sort we’ve just been discussing is V(x) = O every- 
where (the free particle). The even and odd wave functions for energy, E, in 
this potential are cos(kx) and sin(kx) where k = J/2E in atomic units. The 
full solutions, (x,t), for the states of energy, E, are then 


~ik 2/2 ~ik 2t/2 (31) 


(x,t) = cos(kx)e and sin(kxj)e 


We will now put together sets of these states to make pulses, or wave packets. 


FREE WAVE PACKETS 


Since for the continuum we have a largé range of closely-spaced allowed 
energies and wave vectors, we can build up a shaped pulse. Such a pulse is 
usually known as a wave packet. To make a Gaussian-shaped wave packet, 
exp(-(x/x,)?), for example, you use a range of wave vectors, k (around k = 
some central k,), each with its own amplitude, A(k), where A(k)& x, exp 
(-(k-k,)? X97/2). (This is just a formula you can find in books on Fourier 
transforms.) 


You can build up almost any shaped free packet you wish. For a square wave 
packet with width x,, A(k) = X 9 sin((k-k )x_/2)/((k-k, )x,/2). For a trian- 
gular wave packet with basic width 2x5, Atk) = Xo sin? ((k-k 9) X9/2)/((k-kg) 
Xo/2)?. The central wave vector, k,, determines the speed of the packet. The 
group velocity is v = dE/dK = hk,/m or kg in atomic units. 


The next section discusses what happens when you set up such a wave packet 
at time equals zero and then let go. Remember that the wave packet is travel- 
ing through a dispersive medium; the packet shape will distort and deform as 
the higher energy states in the packet shift out ahead of the lower states. 


THE PROPAGATION OF FREE WAVE PACKETS 


A program can easily calculate and add up the amplitudes of the various com- 
ponents of a packet. Because of limitations on how many amplitudes can be 
stored and how long it may take a program to RUN, an approximation to the 
continuous distribution of amplitudes, A(k), must be made. Using twenty k’s, 
for example, across the Gaussian A(k) produces a reasonably good Gaussian 
wave packet in space. 


We set up a packet, centered at the origin, at time = 0; we then let time pass 
and calculate each component's time factor, exp(-iEt). This part of the full 
wave function shifts the relative phases of the various components of the 
wave packet (since each component has its own energy, E). At some later 
time, you add the phase-shifted components back up. The packet has moved 
and distorted. Notice that the time parts of the wave functions do nothing to 
the amplitudes of the states because | exp({—iEt) | = 1. 


Figure 3 shows a block diagram of this strategy. One implementation of this 
strategy is the following program: 


Figure 3. 









Calcutate 
x, Wavefunction 





Flow Chart of Computer Approach to 
Wavepacket Strategy 
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FREPAC 


268 
270 
288 
290 


186 DIM PC21,.26).0C 2). 263, E021),CC 213,50 21) 
28 PRINT “CENTRAL ENERGY, HALF-WIDTH OF PACKET?" 
38 INPUT E8.sX6 
50 LET K9"3/¢106x6)<— Ak 
6@ FOR I2! TO 21<——-—___ 
76 LET K@=KB+(I-11)*K9 Step through states . 
88 LET El13=K@ske/2~<— Energy of each state Store energies 
98 LET K1*K@-K8s 
10@ LET Aw. 3995696X 6eEXP( =x 68x 6eK 14K 172) -—A plitude and wave 
11@ FOR Jai TO 26<——___ Step across X of each state functions for 
128) LET X@=-2.5¢(J-1)8.58 . 
19@ LET PCt.J1=aecos¢(Kexxgo<«—Aeal part of spatial several states 
146 LET QCIs-JIPASSINCKO*#XO) wave function 
156 NEXT J ae 
168 NEXT I Imaginary part 
178 PRINT “TIME THE PACKET HAS PROPAGATED?" 
18@ INPUT TO 
198 FOR 1! TO 21 
208 LET CCIJ=COSCELII*TA) iy 
218 LET Eran: Store exp(-iEt) 
226 NEXT 1 
238 FOR Jeli TO 26 q ~~ Step across X 
248 LET XO=-2.5¢(J-1)#.5 
T 
eae FOR tel 70 21a-—— otep through states 
LET PieP1+K9e (PEI sg) *CC1I-OC1.0)*St1))-«—Real part of y (x,t) 
corre le JI*COlLI+PCIs,JI*SC1I) I~—/ Maginary part 
308 PRINT X@,Pi*Pl+¢P2ePe 
316 NEXT J 
326 PRINT 
33@ GoTo 17¢«—Aeturn for new propagation time 
34@ END 
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GET-FREPAC 
RUN 
FREPAC 


CENTRAL ENERGY, 


72. 1 


TIME THE PACKET HAS PROPAGATED? 


78 
“2-5 
~2 
-1-8 
=-1 
=.§ 
3 
25 
t 
1e5 
2 
2.5 
3 
3-5 
4a 
4.5 
5 
§.5§ 
6 
665 
7 
75 
3 
8-5 
9 
9-5 
16 


TIME THE PACKET 


71 
-2.5 
-2 
-1-5 
=) 
-25 
@ 
5 
I 
168 
2 
- 2e5 
3 
3-5 
a 
4.5 
5 
5.5 
6 
605 


=a 


725 


5 
5 


6 
8 
9 
9 
is 


TIME THE PACKET 


73 
-2-5 
-2 
=1-5 
-1 
-5 
8 
oS 
i 
1e5 
2 
2.5 
3 
3¢5 
4 
4-5 
S 
5-5 
6 
6.5 
7 
7.5 
8 
8-5 
9 
9.5 
16 


TIME THE PACKET HAS PROPAGATED? 
? 


HALF-WIDTH OF PACKET? 


2. 61288E-63 
1-86537E-G2 
«08321 

2 37 88 89 

©78 1667 

le 

2781667 

¢ 376889 
©105321 
1e-80537E-82 
2001 288E-63 
1+ 44626E-84 
te S6646E-86 
7.619 04E~88 
92 5@167E-87 
ie 7588 6E-87 
8 © @367 6E-87 
7.47309 E-68 
7¢27233E-67 
2e83461E-88 
6. 635S4E-67 
7¢7T7TOS3IE-89 
be 15697E-07 
4. 21349E-18 
5.859 46E-87 
16 53338E-89 


HAS PROPAGATED? 


2e 369 3JE<85 
Ze 61658E-04 
7097562E-083 
3e 128S58E-82 
9 $299 3E-82 
«22992 

e 431434 
626177 
«768869 

¢ 626177 

© 431434 

0 22992 
9e52989E-92 
3e 12857 E~62 
7 097 S64E"“83 
1-59562E-83 
2-61 668E-84 
20369 LLE-85 
208 279SE-86 
3¢ 0655 1E-87 
§ 0838 85E-067 
Le 12792E-67 
Se &47S68E-87 
&e 80029 E- 86 


HAS PROPAGATED? 


2- 4689 £E-G4 
4.879 88E-84 
le LIS67E-83 
2e 3168 4E-83 
4.547 1SE-83 
« 088767 
1e53854E-82 
«826273 
Ae28155E-82 
6 37 A96E~-82 
9e SI SB4E-B2 
e 128852 
0 169734 
© 2139225 
2 253389 
¢ 28 6229 
« 389 357 
© 318879 
« 389357 
= 26 6228 
2 253389 
213224 
169734 
e 126852 
9e-I1LIBIE-B2 
6-37 486SE-C2 
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To use terminal plotting you might make the following changes: 


S READ X8.X9,Y6.Y9.Z20 
6 DATA -8eS.s 1626. 1286 

7 GOSUB 9660 

360 LET YO=PiePi+PeePe 
362 GOSUB 9166 

326 GOSUB 9266 
APPEND=-TTY PLO 

RUN 

FREPAC 


CENTRAL ENERGY, HALF-°WIDTH OF PACKET? 
%2. 1 
TIME THE PACKET HAS PROPAGATED? 


6 

MAX Yo 1} 

gg IARI NAAT AAT OANA AMARA NANAK AN MHA INA AA AN ANAAA AANA NR 
g 

Y 

Y 

Y 

Y¥ 

Y 

Y 

Y 8 8 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y @ 6 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 8 8 

Y 

. @ 8 

Y@ ®@8e¢66¢686868080«6G6=~«hUCUmKUCUKC CU CS 


KXXXXKXAXAXAAXAXAXAAXAKAXARRAANAAK KAKRKXAKK XXX AXKX XK XARA KXAN KKM KKKK K KAA KKH 


MIN Y= 8 MIN X=#-2.5$ 


TIME THE PACKET HAS PROPAGATED? 
71 
MAX Y= 1 


MAX Xe 18 


XKXXAXXAXAXKAXAAK XX XAAAAXAXKXAXKX XX KA KKK XXX AAAKK KRAKK KXRXARAKAKA RAK AA X AAKK KAXKAKAK 


Ree eee KKK KK KKK KKK 


@ 
ye @ 289 


MIN Y= 6 MIN X=-2-5 


MAX X® 18 


AXXAAXKRXXKAXXXAXXX 
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TIME THE PACKET HAS PROPAGATED? 

23 

MAX Y* 1 

XNXXKXXKXXAAXAKXAKX AK AX KX XX XXX KKK KAKX KK ARK KXAKKAKARAKK AK KKK KAKKK KAKA KX KK KK KARR 


See KKK KKK KKK KKK 


8 8B 

Ys @ 6@ 8 6 @ 
XXXXXAXXXXXXAXXXKKAAXK MAX RA KAA RA AKKMKK KX AKAKAK AANA AK MK AK X KKK K KK KK RX ARK AKKAKK 
MIN Y= @ MIN X#-2-5 MAX X# 18 


THE PROPAGATION OF A PACKET 
OVER A POTENTIAL WELL 





You can also use this strategy to propagate wave packets by potential wells. 
For a general (symmetric) potential well, you could calculate the spatial wave 
functions numerically as discussed in Chapter One. The strategy of the calcu- 
lation is the same as that for free packets. The only difference is the calculation 


of the spatial wave function. 


We will form our packet at x = 0 to the left of the well; we will then start the 
packet towards the well. In the end, part of the packet will get past the well; 
part will be reflected. Some of you may have seen the films of the quantum 
mechanical scattering of packets off the finite square well. The smoother the 
potential, the easier it is for the packet to get past. Smooth edges on the 
potential also diminish the rapid oscillations (the interference effects) at the 


edges of the well. 
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A program to propagate wave packets by potential wells might look like the 


following: 


GET-VPAC 
LIST 
VPAC 


18 

28 

38 

ae 

cy) 

68 

76 

88 

98 

198 
118 
120 
138 
148 
15¢ 
168 
178 
18 
199 
208 
210 
220 
230 
248 
250 
260 
276 
280 
29¢ 
300 
31 
329 
338 
340 
35 
366 
3708 
382 
390 
Ate 
Al 
420 
430 
Aas 
458 
460 
a7 
486 
a9 
520 
518 
529 
538 
548 
550 
5668 
57 @ 
588 


5908 
{ 689 
6198 
623 


638 
648 
656 
668 
676 


DIM PC21, 262,00 21, 263, EL 213,CL 213.5021) 

DEF FNV(X) =~ 1 GB" EXP( -(X=6) #(X=6) <+———Potential, V{x) 
PRINT “CENTRAL ENERGY. HALF-WIDTH OF PACKET?" 

INPUT E8,X6 
LET K8=SQR( 2eE8) 

LET K9=3/¢ 18%X6) 
Pt aes 9|<«——— Step through states 
LET K@#K8+CIe1l1)*K9 

LET EC 1)=KaeK0/2 «—Energy of each state 

LET K7#K0~-K8 Amplitude of 
LET Ame 399569*X 68 EXP( -X6#X 6#K7#K7 /2) 

LET QCI,1]=Q@=Pia@ each state 
LET PCI. 1}*#PO"A a hads ; 

eer Initialize calcula tion 

LET X0#X529 of radial wave function 

LET N#1 

LET Ki=P] 

LET M1#Q1 

LET LO=2e(FNV(XO)-EL 12) 

LET L1sLa*Pe 

LET Ni=L@*Qe 

LET K2=Pi+L 1lex1/2 

LET MO=QieN 1#X1/2 

LET L@mQeC FNUCXO4X1/2) EC I} 

LET L2=L@u( P@+K 1X 1/2) Fourth Order 

LET N2sLe@*( O6+M1*X1/2) Runge-Kutta 

LET K3@P1+L2eXxX1/2 ; 

LET M3=QL+N2"x1/2 calculation of 

LET L39L0*( P@+K2*X1/2) real and imaginary 


LET N3*LQ"( QG+M2*x 1/2) F radi 
LET K4=P1+L36X1 parts of radial 


LET MAsQloNJex! wave function 
LET LOm2e( FNVCX64X1) EL 1)) 


LET L4#L@*¢ PG@+K 3*X 1) 

LET N&sLQ*(QG+M3*X 1) 

LET P@=@PaexX 1l* CK 14 2eK 2+ OwK 34K 4) /6 

LET Pl#PleXl*(LJ+2eL2¢ 2*L3+L4) /6 

LET Q@2#QG0+X1*(M i+ 2aMee CaM 3+M4a)> /6 

LET Qi@Ql+X1#(N 1+ 2aN 24 26N 34N4) /6 

LET X@=X@+X! 

LET X5=XS+X1 

IF ABS(XS-.5)>X1/2 THEN 488 Store real and 
LET X5a@ g ; 

LET NaNo1 imaginary parts 
LET PLIZNI=PO . . 
cee aed wisKe of spatial wave functions 
coe te THEN 1804 Dorurn for next Ax 

PRINT "TIME THE PACKET HAS PROPAGATED?" 

INPUT TO 
FOR l=i TO 2) 
LET CLIJ#@COSCELII*TO) 
LET SCIJ#SINC-ECII*TO> \ 
NEXT I ) 
FOR Je} TO 26«—__ 

bE KasCd-ij00$ Step across x 

LET Pi=p22g 

FOR anes 2\<——_——Step through states 

LET PIP loK9#( PCI JIeCLII-QC1,d}*SC11)-<«—Fleal (x,t) 

LET P2eP2eKOeCOLI»TI*CEII+PCL,JIeSC1))«q__ 

peers Imaginary (x,t) 
PRINT XO, Pl#PleP2*eP2 

NEXT J 

PRINT ; ; 

goto see~—_—Afeturn for now propagation time 

END 


Store exp(-iEt) 
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RUN 
VPAC 


CENTRAL ENERGY. 
71601 


TIME THE PACKET HAS PROPAGATED? 


76 

8 

05 

‘ 
te5 
2 
2-5 
3 
305 
a 
45 
5 
§.5 
6 
6-5 
7 
7e5 
8 
8.5 
9 
9-5 
16 
18-5 
M1 
1-5 
12 
12.5 


TIME THE PACKET 
71-5 
@ 

5 

1 
1.5 
2 
2.5 
3 
325 
4 
4.5 
5 
5.5 
6 
605 
7 
725 
8 
8-5 
9 
9-5 
18 
16.5 
il 
11-8 
12 
12-5 


HALFeWI DTH OF PACKET? 


le 

° 781693 
0371724 

© 166117 
816311 
2-851 15E-863 
1¢47756E-64 
be 79255E-86 
6.89511E-68 
2e IS27TSE-O7 
3e 11 307E-67 
1e 37568 E-87 
2- 7957 2E-86 
2-55 129 E-89 
Te 4L13IVE-88 
6-6 41582E- 68 
Le 868 62E- 87 
1-98 623E-97 
le 2iSi17E+O87 
1e891854E-86 
2e IDIE6TE~B6 
be 71 39SE-86 
Je 1468 SE~-86 
2e-9TIGTE-G6 
le 72569 E-06 


HAS PROPAGATED? 


20 677 1 3E-86 
3¢3578B1E-85 
$4 418S87E-84 
66599 83E-G4 
2e 39 18 3E-83 
8 - 8447 2E-83 
22 29857E-62 
$6 83398 E-82 
6- 18695E-82 
52 A448 77E-92 
4. 349 22E~82 
3e BL 167E-62 
2- 62889 E-62 
3-847 a2E-82 
7e7TO667E-82 
e 141367 

e 17333 

¢ 168835 

© 121496 
7+«66528E-62 
308 S898 E-62 
1-65678E-62 
5.8877 1E-83 
10717 66E-83 
42 S9 A41E-84 
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Using terminal plotting you get: 


CENTRAL ENERGYs HALF-WIDTH OF PACKET? 
718, 1 
TIME THE PACKET HAS PROPAGATED? 
29 
MAX Y= } 
KNIKMHMHK NHK KK HK HIN KK HNN NHK M KH MMH NH MM MK IMM KH HH MH HK H RK HHH KHHK HK K MMMM HK KHHK KAA HK KKK 
yy! 


I 


SKK SKK KKK KKK KKK KKK KKK 


!o@aeka&: £€ tao 4 do a toy to § 8 dd 
AXAXXXARKAKAARAKX AX AKA RAXAKAAKAAAA NK NAKA RN NK RAHA RAHA KK AAK RK KH AKA HK KKH HAA NAKA AA AN 
MIN Yo 6 MIN Xe @ MAX Xe 12.5 


TIME THE PACKET HAS PROPAGATED? 

T1eS 

MAX Y= | 
XNXAXXXXXAAXAAAKAKAAKR AAR KAN AANAMKAANAAANKKA KKK MAA AA AAAA AA ANAK A MMA AAKMMAK 


1 tft 1 1 | 

t ! t1 1 ] 

Yi adedtsh?:d 41 tt 41 3 
MXXXXAXXKKKXAXAKKRAN RAK KX XXNAK X KANN NHK XX KKK XK MAK KK XK KK KK KK KX KXAN AK KKK KA KAA 

MIN Y= @ MIN X= @ MAX X= 12-5 

TIME THE PACKET HAS PROPAGATED? 

? 

END 


KRM eK RK KKK KKK KKK KKK KKM 


The program uses a fourth-order Runge-Kutta method of integration to find 
the spatial wave function. Adding up wave functions to make packets demands 
that not only the amplitudes of the wave functions but also the phases must be 
very accurate. You must be careful to use a small enough space step. 


Because of the limitations of some computers, only 21 different wave vectors 
are used in the program presented. Sums over only 21 wave vectors are hardly 
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integrals over all wave vectors. This compromise leads to ‘‘ghost’’ packets 
which image the main packet but lie at large positive and negative x. The com- 
promise also means that the packet shape is only approximate. This last pro- 
gram takes along time to RUN. On some machines it is a half hour before the 
program finishes finding the spatial wave functions. 


EXERCISES — Use computer when appropriate 


!. Using the fact (from Fourier Series calculations) that a triangle wave has 
only odd harmonics in the ratios of 1/{2n+1)? n=1, 2, 3, ...: 


7. Set up such an initial electron distribution in an infinite square well. 


2, Find the shape of the combination of states after 1 through 5 units 
of time. 


3. Identify times when certain modes seem dominant. Explain this 
phenomenon in terms of the periods of the modes you used. 


"Il. Using the fact (from Fourier Series calculations) that a linear periodic 
wave (a ramp or sawtooth wave) has harmonic content in the ratios +1, -1/2, 
+1/3, -1/4, ... : 


7. Set up such an initial electron distribution in an infinite square well. 


2, Find the shape of the combination of states after 1 through 5 units 
of time. 


3. Identify times when certain modes seem dominant. Explain this 
phenomenon in terms of the periods of the modes you used. 


“ltl. Form a combination of the 3 lowest states of the infinite square well. 


1. Find the form of your combination for T = 0, 7, 2,9, and 16 units 
of time. 


2. Is the motion periodic in time? What is (half) the lowest common 
multiple of the periods of the states? 


IV. Form a Gaussian free wave packet at the origin at time equals zero. Use a 
central energy of 8 and a packet half-width of .5. See figure at right. 


1. Propagate the wave packet for 0,.5, 1, 1.5, 2, and 2.5 units of time. 
Observe and discuss the behavior of the packet, 


2. From your calculations, measure the group velocity of the wave 
packet. Compare this result to the central wave vector of the pulse and 
to the phase velocities for the lowest, highest, and central energy com- 
ponents of the packet. 





3. From your calculations, measure the full-width at half maximum for 
the packet at each time. The width is predicted (for a Gaussian packet) 
to go as W(t) = ((W(0))2+ t2/(W(0))2). Do your measurements of the 
propagated packets agree? 


4. By integrating (crudely) ||? across the whole range of x, see if the 
total probability for the particle remains constant. 


“V. Form a triangular free wave packet at the origin at time equals zero. Use 
a central energy of 8 and a packet half-width of .5. See figure at right. 


1. Propagate the wave packet for O,.5, 1, 1.5, 2, and 2.5 units of time. 
Observe and discuss the behavior of the packet. 
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2. From your calculations, measure the group velocity of the wave 
packet. Compare this result to the central wave vector of the pulse and 
to the phase velocities for the lowest, highest, and central energy com- 
ponents of the packet. 


3. From your calculations, measure the full-width at half maximum for 
the packet at each time. The width is predicted (for a Gaussian packet} 
to go as W(t) = (W(0))2 + t2/(W(0)2)”. Do your measurements of the 
propagated packets agree? 


4. By integrating (crudely) | ¥ |? across the whole range of x, see if the 
total probability for the particle remains constant. 


*VI. Form a square free wave packet at the origin at time equals zero. Use a 
central energy of 8 and a packet half-width of .5. 


1. Propagate the wave packet for O,.5, 1, 1.5, 2, and 2.5 units of time. 
Observe and discuss the behavior of the packet. 


2. From your calculations, measure the group velocity of the wave 
packet. Compare this result to the central wave vector of the pulse and 
to the phase velocities for the lowest, highest, and central energy com- 
ponents of the packet. 


3. From your calculations, measure the full-width at half maximum for 
the packet at each time. The width is predicted (for a Gaussian packet) 
to go as W(t) = ((W(O))? + t?/(W(0)? }7. Do your measurements of the 
propagated packets agree? 


4, By integrating (crudely) | W |? across the whole range of x, see if the 
total proability for the particle remains constant. 


**VII. Usinga Gaussian wave packet (of central energy E = 50) and a Gaussian 
well (of depth 100 and half-width 1): 


1. Propagate the packet for times of 0,.2,.4,.6,.8, and 7. 


2. Describe and discuss the results in terms of reflection coefficients, 
transmission coefficients, and interference effects. 


**VIIL. Using a triangular wave packet (of central energy E = 10) and a 
Gaussian well {of depth 100 and half-width 1): 


1. Propagate the packet for times of 0, 1, 2, 2.5, 3, 4, and 5. 


2. Describe and discuss the results in terms of reflection coefficients, 
transmission coefficients, and interference effects. 
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CHAPTER FOUR: THREE-DIMENSIONAL CONTINUUM STATES 
- FOR SPHERICALLY SYMMETRIC POTENTIALS 


INTRODUCTION 


With the computer you can solve for the continuum states of any spherically 
symmetric potential. The method is very similar to that described for caicu- 
lating the discrete states for these potentials (Chapter Two). The three- 
dimensional time independent Schrodinger equation is 

h2 = 

~37 Veer) + VT) YT) = Ey(T) (32) 

m 
As we did in Chapter Two, you can separate this differential equation into 
polar coordinates when V is spherically symmetric, i.e., when V(r) = V(r). 
You set W(r) = R(r)Y¥(6,¢) and substitute this expression into Equation (32). 
The result is (1) an equation for R(r) which includes the effective potential 
and the energy, and (2) an equation for Y(@,@) are spherical harmonics: 


VEO 0) Noe (cos@) e=!M? (33) 


where No is a normalization coefficient, 2 and m are integers such that 
l\m|< =, and Lo is an associated Legendre function. 


The equation for R(r) can be simplified by transforming to an equivalent for 
P(r) =rR(r). The equation for P(r) looks very much like the one-dimensional 
Schrodinger equation: 


d2P(r) 


sa = AW Ir) + QL 1) /(2r2) ~ E) Pir (34) 


This equation can be solved iteratively in the manner discussed in Chapters 
One and Two. 


BOUNDARY CONDITIONS FOR P(r) 


P(r) must go to zero smoothly as r goes to zero. This is just a consequence of 
the fact that R(r) must be bounded everywhere. tf R(r) were infinite at some 
fy, then |v (r, 6,9) [2 would also be infinite at that point. This would mean in- 
finite probability of finding the particle at tr, which is unphysical, 


This demand that P(r) go to zero as r goes to zero demands that the solution 
for P(r) be unique. Hence, unlike the case of continuum states in one dimen- 
sion (Chapter Three), there is only one radial wavefunction for each continuum 
energy in three dimensions. 


If V(r) can be approximated by some simpler potential (for example, the 
Coulomb potential for the hydrogen atom) for small r, then you can approxi- 
mate the behavior of R(r) in this region by that of the solutions to the simpler 
potential. Then the behavior of P(r) is also known for small r. (For the Cou- 
lomb case P(r) ~ r£*1), This allows you to initialize the wave functions for 
such cases near r = 0 and then integrate step-by-step out to large r. When this 
is possible, it is the easiest strategy to use. 


Assume now that V(r) goes to zero for large r. Most important, three- 
dimensional potentials do go to zero far enough from the force center. Then, 
as you can see from Equation (34), the form of the Schrodinger equation for 
P(r) becomes P’’ = -2EP. The continuum states for such a potential have posi- 
tive energies, E, so this equation for P has the (asymptotic) solution P(r) = 
Asin (kr) + Bcos(kr) for larger, k = JE. Since there is only one state for each 
energy, we can conveniently write this solution as P(r) = A sin(kr-27/2+6 9) 
for large r. dg is called the phase shift. The equation for R(r) at large r is then 
spherical Bessel functions, R(r) = sin(kr-&/2+6o) /r. 


52 UO QUANTUM MECHANICS 


Notice there is no reason that P(r) should go to zero as r gets large. This is un- 
like the situation for discrete states; in discrete states, the wave function 
entered a classically forbidden region for large r, and P(r) had to go smoothly 
to zero. For continuum states, the only demand on R(r) and P(r) for large r 
is that they remain bounded. 


THE ITERATIVE PROCEDURE 


You can solve the equation for P(r) (Equation (34)) iteratively in much the 
same way we solved the same equation for discrete states in Chapter Two. 


Suppose you know P(r) at some point r, and the behavior of P(r) somewhere 
else (say at r). You can initialize P(r) and P(r) at r,; you can then take a 
small step Ar towards r and find the next values of P and P’ from 


P{rtAr) = P(r) + P(e A)ar (35) 
ts, } w Ar 
P'(rtAr) = P(r} +P (r+ Sar (36) 


The value of P” is obtained from Equation (34). The use of the half-step 
(r+Ar/2) is discussed in the Appendix. 


In this manner, you walk your way in (or out) r until you reach r,. If the be- 
havior at r, is correct, you have the state for that energy. 


If you do not know the asymptotic expansion for P(r) near r = O (that is, if 
V(r) is not approximated by a potential whose solutions you know), then 
initialize P(r) to sin(kr-&2/2+6,) at large r; integrate back to r = 0; and guess 
phase shifts, 50, until you find the one for which P(r) goes smoothly to zero. 
One possible implementation of this strategy is: 


SCATI 


16 DEF FNVCR)#46@e(Rt(-12)-Rt 6-6) -«——_ Potential, V(r) 

26 LET Ps3.14159 

36 PRINT “RCSTART), RC END)» DELTA-R?" 

4@ INPUT R9,s RB, R7 

$@ PRINT "ORBITAL #(L). ENERGY, PHASE DEGS. ) ?" 

68 INPUT L,E,D 

7@ LET DepeD/186 a 

88 LET K#SQR( 2eE)«——— Wave vector tializati 
SG: ea eee cane cay, [Initialization 
100 LET P@=SINC(CK*RB+D> 

110 LET PleK*COS(K*R8+D) 

1286 LET R=R9 

13@ LET N=e 

14@ PRINT R, P@,P6/R 

158 LET RS@R9+(RE-RI) /26 


168 LET R6=R+R7 /2 e 
176 LET P2e2eC FNVCR6) +Le( Le 1) /( 2e@R6eR6)-E) ®& (P84 P1*R7/2) Int grate 
19@ LET P3=P6 Schrodinger 
196 LET P@=Pd+(P1+P2eR7 /2)*R7 : 

208 LET Pl=Pi+P2eR7 equation 


2t@ LET R=R¢R? 
228 IF SGNCPQ)*SGN(P3) THEN 269-<-———— /est for zero crossing 
239 LET NeN+1<———- Count zero crossings 

248 IF N>1 THEN 268 \ 

250 LET R@=((R-R7)*Po-ReP3) /cPa-P3> f Otore asymptotic 

268 IF R>RS THEN 294 j zero crossing 

279 LET R5=R5+(RB-R9)/28 Print loop 

260 PRINT R,POsPO/R 

29@ IF R>R8 THEN 169-<————Aeturn for next Ar 

300 LET P9=N«P-KeRO+LeP/2 \ , 

31@ PRINT “PHASE SHIFT (DEGS.)> ="s1aeep9/pf Print phase shift, 5¢ 
328 PRINT 

33@ GOTO 3¢<—Feturn for new parameters 

348 END 
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The parts of the program dealing with phase shift calculations will be discussed 
in a section below on Scattering Theory. 


If you do know the small r behavior of P(r), then the solution is completely 
determined. You initialize the wave function P(r) near r = 0 and integrate out 
to large r. You can then measure the phase shift, 6, , from the positions of 
zeros of P at large r. At a zero, Fg, in the large r region, krg-&7/2 + 5) = n7, 
where nis the number of zeros of P(r) between r = QO and r = fy; k = E. rg is 
measured so 6, can be calculated. n7 is the actual phase of the wave function 
at that zero; kry ~ &7/2 is the phase the wave function would have if V(r) were 
zero everywhere (i.e., when P(r) = sin(kr-R7/2) everywhere); so 5, is the 
phase shift of the true wave function relative to a free particle. Let’s look 
at a specific problem—continuous states for the screened Coulomb potential, 
V(r) = -exp(-r/8)/r. 


CONTINUUM STATES FOR A SCREENED COULOMB 
POTENTIAL 


You already know (before you start) a good deal about the solution. Since the 
potential is approximated by the unscreened Coulomb potential for small r, 
the wave function, P, comes out of r = Oas r2*" ; the derivative then behaves 
as (2Q+1)r®. At large r, P(r) becomes a sine wave of constant amplitude. tn this 
sine wave region, we will measure the phase shifts for a given energy and vari- 
ous orbital angular momentum quantum numbers, &. 


Find the wave function of the state with E = 4 and & = 0 using the computer. 


A program implementing the strategy we just discussed follows: 


SCATe 


1@ DEF FNV(R)=-EXP(-R/8) /R-«— Potential, V(r} 

28 LET P#3.14159 

36 PRINT "RCSTART)sRCEND), DELTA-R, ORBITAL #@(L). ENERGY?" 

4@ INPUT RB,RO-RTI-LVE 

$@ LET K#®SQR( 2*E)-«—___- |W/ave vector 

68 LET PO=RBI(L41) ee 
7@ LET Pl#(Lo1)*RBtL Initialization 
8@ LET R=RS 

98 PRINT FR. P@.PO/R 

108 LET Nag 

11@ LET RS#CRO-RB) 720 | 

120 LET R6sR+R7T/2 

13@ LET Peseec FNVC RG) +L *(L4 1) /C 2 R68R6) -E) *(PO+P1*R7/2) Integrate 


146 LET P3=Pe | rodin 
15@ LET P@sPB+(P1+P2eR7/2)*R7 Sch o gor 
168 LET Pi=plePp2eR? equation 
178 LET R= ReR7 


188 IF SGNCPO)=SGNCP3) THEN 21¢<—Jest for zero crossing 

19@ LET N=Ne! 

208 LET R@=((R-R7)*PO-ReP3) /(PB-P3) Count zero crossings 

21@ IF R¢RS THEN 248 | p.. Save fast zero crossing 

226 LET RS=R5+(R9-RB) 728 | Print loop 

238 PRINT RPG, PO/R 

249 LF R<R9 THEN 120<———Aeturn for new Ar step 

258 LET POsSNeP-KeRO+LeP/2 ‘ : 

26@ PRINT “PHASE SHIFT (DEGS.) ="s 186*P9/P | Print phase 

ae PRINT | shift, do 
GOTO 38 

298 END “—— Return for new parameters 
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SCAT2 
RC START). RC END). DELTA-R,- ORBITAL #(L),. ENERGY? 
PTIE@$82 204-8150. 4 


PHASE SHIFT (DEGS-) 2 76.8723 
RC START)» RC END)» DELTA-R, ORBI TAL 


1. @BOBGE- 88 1» G@SGQRE-28 1 
101 -. 104641 ~. 193685 
2.81 8. 31G46E-02 4. 13456E-82 es 
3.20802 -3.26691E-82  -1.08896E-G2 
4-98804 -2022426E-82 -5-S6858E-83 
5.08887 8.83536E-02 1. 68785E-82 
eT) ~.135873 ~2.25118E-82 
7.08013 + 186357 2.57649 E-22 
8.60014 ~.21112 ~2.63895E-@2 
9.00816 + 223789 2.48 561E-G2 
18.0882 =. 216279 -2.16275E-@2 
11-0862 - 188997 1.718 12E-92 
12.8802 ~. 144843 -~1.20033E-82 
13.9992 - 085408 65569 72E-83 
14.0903 -1685095E-02 -1.322@8E-023 
15.9803 ~5e03436E-82 -3-35618E-e3 
16.8863 0114574 7. 16673E-83 
17-8001 =~. 16799 -9 688 166E-83 
13-81 2027 69 1.1258 7E-92 
19.6998 ~. 222463 ~ 1. 16994E-B2 
28-8896 «22866 1. 18277E-92 


PIE"“8, 2052612144 


PHASE SHIFT (DEGSe) = 58.1798 
RC START), RCEND) » DEL TA-R,» ORBITAL 


#(L), ENERGY ? 


1+ O@8BBE- 08 1+ BBAGBE-16 1+ GBBBAGE-O8 
1-81 9+9S5384E-67 9-85450E-07 
2-65 ~1e89195E-06 -5¢43257E-87 
3+ 88862 $6 11603E-06 3e 72086E-87 
4-8E284 ~1-85186E-06 ~2e62762E-87 
S- 88887 8+99 266E-07 1-79851E-67 
6+ 88899 ~6+¢67755E-87 “121129 1E-67 
7+886811 3-73776E-67 5-339 $8 E- 88 
8+ 88814 ~4- 23069 E- 68 -5-31327E- 89 
9- B8816 ~2-958 26E-87 ~J- 28 69 2E- 88 
18-@882 6. 89 248 E-07 60 B9237E-88 
11-8882 ~8 - 67 447 E-87 -7-8857 4E- 28 
12-8802 le C4499 E-86 8-788 @SE-88 
13-8092 ~le 12398 E-86 -8- 6458 1E-88 
14-8883 1-89 68GE-26 7+8 28 4BE-08 
15¢8883 -9« 63826E-87 = 6¢ 42684E~- 88 
16-0883 737 329E~-87 4- 688 22E- 88 
$17 - 8881 ~4.48326E-087 -2eS981I3E- 68 
18-981 $e 32242E-07 Te 34273 E- 89 
19-8098 2@-IBLILE-87 le LA736E-88 
20-9896 -5-A7661E-07 -2-73668 E- 88 


PIE"8s 28s 681s 254 


#(L)s ENERGY? 


1- 82008E-88 1+ @BGAGGE-24 le BBCSGE~ 16 
1e@1 4e58785E-12 4e-S4243E-=12 
2-41 -~1+-86842E-12 “92 29S564E=-13 
3088862 ~1¢694727E-13 -3- 49 889E=14 
4-88804 1. 69441E-12 4e 23598 E-13 
5-898087 ~3-83941E-12 ~6¢@67874E-13 
6- 88289 4eB6168E-12 6.769 24E-13 
7+88011 ~4. 68843E-12 -6+ 63622E-13 
8-88814 4-842@6E-12 6¢O35247E-13 
9- 88816 ~“4eS3217E=-12 ~5-@356SE-13 
18-8082 3e78OS9E-12 3¢ 78BS3E=13 
Lie 8882 -2- 65989E<12 “20-418 63E-13 
12-6882 le 2788 6E-12 1-86578E-13 
13+ 8982 2° 27821E-13 le 75243E-14 
1426883 “1.7 1282E-12 ~f<e22342E-13 
15.8993 3¢B3I@6SE-12 2. 8204BE-13 
16-6883 ~4.65265E-12 -26S3248E=-13 
17-8881 4- 67676E-12 2-7TSIGIE=13 
18-81 -4-8 479 4E-12 ~2-69181E-!3 
19-8398 4-58227E- 12 ee ALB4TE-1t3 
20.8896 ~3e86771E-12 -1.93292E-13 


PHASE SHIFT (DEGS-) = 58.2172 
RC START)». RC END). DELTA-Rs ORBITAL #¢(L),. ENERGY ? 
PlE-8, 20,-815304 


1. 8BSEBE- 38 le BABSBE-32 1. @@BBGE-24 
1-8 9-29985E-18 9+ 20776 E- 18 
2-81 1e39632E-17 6-9 468 6E=- 18 
3- 88882 ~I1e87I@1E-17 -6. 24333E- 18 
4-80664 1e85@72E-17 Ae 6267 4E- 18 
5- 88887 -~1-S7T1S8E-17 ~3- 14311E-18 
6 68889 le 11726E-17 1-86208E-18 
7e@6O11 “-S-SI514E-18 -7-878 64E-19 
8-08814 ~6+ 24836E~ 19 -7-¢88032E-28 
9288916 6- 61848E- 18 Te 3447 6E-19 
18.8802 ~1. 18 457E-17 ~1.-18455E-18 
13-8802 1¢$8172E-17 1+ 4379 8E-18 
12-8882 “1e81431E-17 ~“leSLISGE-18 
13-8882 1-8 608 1E-17 1e43137E-18 
14.6883 -1-71823E-17 -1.-22728E-18 
15-8863 le 48227E~-17 9+ 348 29 E~ 19 
16-8883 9. A5674E-18 “509 1B34E=-19 
17-8881 32948 49 E-18 2+ 32L62E- 19 
18-9! le 4227 6E- 18 7e89985E-28 
19-8898 “7-58 !16E-18 -3. 76S51E-19 - 
20+ 8896 1e21736E-17 6e BB 38 7E-19 


PHASE SHIFT (DEGS-) = &6.7426 

R¢( START). RC END)» DELTA-R, ORBITAL #(L).» ENERGY ? 
? 

END 
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Aterminal plot of this wave function can be obtained by using the subroutines 
in TTYPLO (as discussed in the Appendix). The changes you might make and 
a RUN are: 


S READ X8,X9sY8sY9 
6 DATA 6, 28,-215, 1 
7 GOSUB 9808 

8 LET X9=Y@"Z 020 

9 GOSUB 9186 

99 LET Z@=1 

91 LET X@=R 

92 LET Y9=P8/R 

93 GOSUB 9186 

230 LET X@#R 

231 LET YO=PO/sR 
232 GOSUB 9188 

278 GOSUB 9288 
APPEN D- TTY PLO 


228 LET R5#R5+(R9-RB)/78 
RUN 
SCAT2 


RCSTART) s RC END)» DELTA-Rs ORBITAL #@(L)»s ENERGY ? 

7 1LE<“8s 2829315024 

PHASE SHIFT CDEGSe) ® 76-8723 

MAX Y= | 
MXXKKXXXKXXKKXXKAKKAKKKXKAA KKK KK AK KK AK KK KK AHN K KKK KKK KKM K KK NMA HK HKHK KX NHK KAKK KH 


<< 


I Hl 11 i 
D2 ER PPR ADD TEU DDEDELDDEDDA DDD 


Ba a aaa a a i a en a i 


XXKAXAKXKAXAXKAAKAKAAAKK AA AKA AK KAKAKAAKAK AXAAK KA KRXAKAK AK K HK A ANAK KK XK AA KAA XH KK KK 
MIN Yeo. 15 MIN X= 6 MAX X® 286 

R¢ START), RC END), DELTA-R, ORBITAL #@(L)» ENERGY ? 

? 


END 


SCATTERING THEORY 


The presence of a potential V(r) modifies the wave function (compared to 
the wave function for the free particle). Scattering theory deals with the 
phase shifts tn terms of the differential and total cross sections. 0(@) and 
Oro¢ - Using the ideas developed in the last section, you can use the computer 
to find these phase shifts and cross sections for general, spherically sym- 
metric potentials. 


First, (to make clear the meaning of the phase shift) let’s RUN the program 
for acontinuum state of the free particle (V = 0 everywhere). You are actually 
producing the free particle wave function as expanded in polar coordinates. 
The radial function R(r) will be a spherical Bessel function. 


We can use the program which initializes the wave function and its derivative 
at small r. For 2 = O the differential equation for P is P’’ = -2EP which gives 
the only acceptable solution, P = A sin(kr) where k = /2E. This solution is 
true everywhere so we already know that R(r) = A sin{kr)/r for & = 0 states. 
Notice that the small r behavior of P is P~ r&*1 . 
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For & >= 1, the Schrodinger equation for P becomes for small r 
PP(r) = (2(R+1) /r2) Pir) (37) 


The solution for P has P ~ r2+?, So these solutions can be initialized at small r 
to P ~ 241 P’ ~ (2+ 1)r®. (These are the same arguments necessary to justify 
this initialization for the Coulomb potential; these arguments work for any 
potential which is dominated by the angular momentum term, 2(&+1)/(2r2), 
at small r.) 


A RUN of the program for V = 0, E = 4, and 2 = Ois: 





GET-SCAT2 

16 DEF FNVCR)«¢ 

RUN 

SCAT2 

RC START)» RC END), DELTA-R, ORBITAL #(L).s ENERGY? 

TIE“8, 2802 -B1,024 

1-O8RB0E-88 1-2@0800E-28 i 

1281 9-IS2ISE-82 908338 4E-22 
2-81 ~« 198986 -2 898998 

Je 28882 « 265323 Fe SIBGE-82 
4c 88884 ~2 3957 38 “82 I93I3I5SE-82 
5.28867  ¢ 383513 7eB7B17E-82 
6- 88889 ~2 336922 ~S-. 61528 E-@e2 
72900811 © 287577 4-18817E-82 
8.88614 -22198275 2. 628 I9E-82 
Fe BBO16 #11283 le 25831E-G2 
18-2882 = 3-8 459 4E-83 #3.84587E-84 
11-8882 -.- {85212 9 e 56454E-63 
12-8882 * 2804842 le 7883 1E-82 
13e 8802 -« 2863635 ~2e 1771SE-82 
14.8893 e 394514 2e 389 34E-82 
15-6983 -2 353472 ©2235643E-d2 
16 8883 © 338867 2e 1128S E-82 
17-8@@1 0 289 797 -1+ 70467E-82 
18-81 221241 le 228 43E~Gg2 
19-8898 -« 128566 60 685 33E-83 
28-8896 1e76864S5SE-82 §-83797E-84 


PHASE SHIFT (DEGS-) = «164799 


RC START) » RC END). DEL TA-R, ORBI TAL 


? 
END 


A terminal plot of this state is: 


16 DEF FNVC R28 


RUN 
SCaTe2 


#(L) » ENERGY 7 


RC START)» RC END)» DELTA-R. ORBITAL #(L)s ENERGY? 

2? 1E-8, 20,281,024 

PHASE SHIFT (CDEGS-) = « 164799 

MAX Y= 1 - 
XXXAAKARAAKAAAXAKAAKAARA KKK KAKA AK AAAAAAANA AKIN KR ARAAKARAARARMMMMM MAM AMA MM MM 
¥i 


| 111 11 11 i i 
@ 4 I A f Pook Gh Fd 1d 2d EAD ELD PETTTADELTAIY 
i Lit Mt i] ti j 


SKK KKM KKK KKK KKK KKK KKK 


XXX XXXAXXXXXXAXAKAKXAKKA ARK KK AXKAAKAXXARKRKANR KA AAAKA AAA KARA AAA MHA KAKA A RXAK KKK RAK 
MIN Ys#-.25 MIN X# @ MAX X= 26 

R¢ START)» RC END)» DELTA-R» ORBITAL #(L).» ENERGY? 

? 

END 
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‘You can see the relationship between the V = 0, & = 0 state and the 2 = 0 
state shown before for the screened Coulomb potential. First, the behavior 
of R{r) at small r is very s-state-like for the screened Coulomb potential; 
the small r region is very reminiscent of the discrete spectrum s states for 
this potential. Second, there is a clear phase shift between the free particle 
and the screened Coulomb wave functions at large r. This difference in phase 
is the phase shift, 6), so central to scattering theory. 


You can calculate 6, directly from your results: Count the number, n, of 
zero crossings of the wave function between the origin and some zero cross- 
ing at rg in the asymptotic region. The phase shift, 59, is defined in scatter- 
ing theory by setting the asymptotic behavior of P(r) = A sin(kr-&7/2+6 9), 
where k = «/2E. (kr-&2/2) would be the phase of the free particle at rg SO 
do = nt-kr, + 7/2. 


Scattering theory points out that you can measure the shape of a potential 
by observing the interference pattern between a particle scattered off the 
potential and free particles which were not scattered. The parameters which 
are important in the interference pattern are the phase shifts, 69, defined above. 


In scattering theory, you write the asymptotic wave function as f(@) exp(ikr) /r. 


The differential scattering cross section, o(@), is then o(@) = |f(@)|2. By ex- 
panding everything in spherical harmonics you find that 


£(8) = 


x|— 


a (22+1) exp (ido) sind, £o(cos@) 
2=0 


where 59 is the phase shift and Lp(cos @) is the &th Legendre polynomial; k is 
still V2E. 


You can calculate the Legendre polynomials most easily by using the recursion 
relation 


RL£o(cos#@) = (2k-1) cos@ L,_, (cos) -(&- 1) L£-5(cos8) 
with the two lowest order polynomials, £, = 1 and £, =cos @. 


A program to do this calculation (and its RUN) might look like: 


18 PRINT “ANGLE (DEGSe) ?" 
20 INPUT T 

30 LET 1T#3.14159eT/18@ 

4@ LET c=CcOoscT) 

5@ LET Loe | <——_—_—_ Lo 

68 PRINT * 6",Le 

7@ LET Lis¢ +—______ 

8@ PRINT " 1%L1 1 

98 FOR L=2 TO 10 

10@ LET Lee c( 24L-1)#CeL!-(L-1)*L@) /L <«— i ; 

ae. Ue eects Recursion relation for L£, 
120 LET L@eL} 

13@ LET LieLe 


146 NEXT L 
158 PRINT 
16@ GOTO 14@~<«——- Return for new angle 
178 END 
END 
RUN 
LEGEN 
ANGLE (DEGS.)? 
745 
@ t 
1 - 7871087 
2 250881 
3 -+ 176776 
a - 248625 
5 ™e 375651 
6 ~+ 148.439 
7 - 127857 
8 29834 
9 + 285536 
16 e11Si14 
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You can measure the phase shifts, 6,, for a given potential and given energy 
from your computer-based calculation of P(r). You can calculate the Legendre 
polynomials using the program given earlier. Finally, you can calculate f(@) 
and, from f(@), o(@), the differential cross section. 
You can calculate the total cross section, 0,,,, even more easily. The total 
cross section (from scattering theory) is given by 


to 


Oror = 4a D7 (2e+1) sin25, 
k* Q=9 


Orr 1S very easily measured from the computer derived phase shifts. 


To calculate cross sections you need to calculate phase shifts for a number of 
orbital angular momentum quantum numbers, &. You need phase shifts for all 
L up to that for which the phase shift is negligibly small. Scattering theory 
suggests that, if R is a characteristic radius of the potential, then: when 
kR >> 1, 69's for ail £ are needed; when kR << 1, only 5, and 6, are impor- 
tant, k =./2E. This fact makes these calculations well suited to individual pro- 
ject work. You could also divide up the different & states and have some 
people find certain phase shifts and other people find others. Everyone can 
then use all the results to derive cross section information. 


The simplest way to calculate phase shifts using the continum wave function 
program is to remove the lines which print out the wave functions. Since 
phase is very sensitive to small inaccuracies in the calculation, the phase 
shifts are only good to a few degrees. You can check the method by calcu- 
lating phase shifts for V = 0 everywhere. You could use a higher convergence 
method for-more accuracy. You would still need to integrate far out to get 
the phase of the asymptotic wave function correct. Even for V = Q (excluding 
2 = 0), you must integrate far enough out to make &(2+1)/(2r2) completely 
negligible with respect to the energy. Notice that, if V(r) becomes negligible 
then both the free particle and the particle in V(r) are obeying P’’ = 2(2(&+1)/ 
(2r2)-E)P. You can get an approximation to the phase shift by comparing the 
phases of the V = 0 and V(r) particles when V(r) is negligible (even though 
neither solution is truly in its asymptotic region). 


For example, find the & = 0, 1, 2, and 3 phase shifts for the screened Coulomb 
potential V(r) = -exp(-r/8)/r. Use an energy, E, of 4. 


The results of RUNing the program are: 


59 


57.1; 


76.8- 0.1 = 76.7; 5, = 58.2- 1.1 


It 
I 
II 
II 


55 50.2 ~- 3.1 47.1; 6. 46.7- 6.4 = 40.3 

If we take the screening parameter, 8, as the measure of the size, R, of the 
potential, then, for E = 4, kR = 22.6 which is farge with respect to &. We 
expect to need a large number of & to produce a good approximation to the 


cross sections, Other kR are illustrated in the exercises. 


Now let's calculate an approximation to the total! cross section using the 4 
phase shifts calculated above. 
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A program to combine the phase shifts into cross sections might look like: 


CROSS 


19 PRINT “MAX. L. ENERGY, ANGLE (DEGS.) ?” 
20 INPUT NDELA 

30 LET A™Ax3.14159/16@ 

46 LET K=#SQR( 2*E) <«—-— [/ave vector 

S@ LET C=COSCA) 

6@ LET St#FisF2=g Initialization 
78 PRINT “DELTA FOR EACH L (DEGS-)?™ 

8@ FOR L=@ TON 

96 PRINT "Le"3L 

16@ INPUT D 

110 LET D=D*3. 141597182 

120 IF Le@ THEN 182 

13@ IF Lel THEN 20¢ 
148 LET L2e¢(2eL~1)*CeL i=(L-1)4L@) A «— L, by recursion 
150 LET L@aLl 

168 LET LisLe 

17@ GOTO 218 

182 LET LoeL2u1<«— L, 

i98 GOTO 218 

200 LET LisL2ec~«+— Jf, 

21@ LET S=SINCD) 

22@ LET FIsF1le¢ 2eL+1)*COS(D)*S#L2 a Freal part of (0) 
238 LET F2eF2+( 2ebe1)«S#S#L2 «-—-———____ /maginary part 
248 LET S1®S14+¢€ 2ebL41>#Se5 — ; 
258 NEXT L Total cross section 
260 PRINT "FCTHETA)=” 

270 PRINT FI/K,F2/K 

286 PRINT "SIGMACTHETAD=" 

298 PRINT (FI*F1+F2eF2) /(K*K) 

306 PRINT “SIGMA TOTAL=" 

330 PRINT 4*#3-.14159"S1/(K*K) 

328 PRINT 

33e@ GoTo 18e<— Aeturn for new parameters 

348 END 


Finally, let’s calculate an approximation to the differential cross section, o(@) 
at 6 = 0 degrees. 


MAX. Ls ENERGY, ANGLE (DEGS.)? 
73. 4,8 
DELTA FOR EACH L ¢(DEGS-)? 
L= 8 
27 6e7 
Ls ] 
787.1 
L= 2 
7aTe} 
Ls 3 
746-3 
FC THETA) = 
2 66521 32 86652 
SIGMAC THETA) = 
16.5969 
SIGMA TOTAL#= 
13+«6242 


MAX+ Ls ENERGYs ANGLE (DEGS.)? 
? 


END 


EXERCISES — Use computer when appropriate. 
|. For the potential, V(r) = --exp(-r/8)/r, and an energy of E = 10: 


1. Plot the effective potential for 2 = 0, 1, 2, and 3. What is kR for 
this situation? 


2. Plot the radial wave function for 2 =O and 2 = 7. 
3. Find the phase shifts for 2 = 0, 1, 2, and 3. 


4, Calculate an approximation to the differential cross section at O, 90, 
and 180 degrees using the four phase shifts from Part 3. 


5. Calculate an approximation to the total cross section using these 
four phase shifts. 


Hint: Use the program which initializes the wave function at small r. To 
get better accuracy for the phase shifts, subtract the corresponding 
values of the phase shifts for V = 0. 
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Il. For the potential, V(r) = -exp(-r/8)/r, and an energy of E = 1: 


1. Plot the effective potential for 2°=0, 1, 2, and 3. What is KR for this 
situation? 


2. Plot the radial wave function for 2 = Oand & = 7. 
3. Find the phase shifts for = 0, 1, 2, and 3. 


4, Calculate an approximation to the differential cross section at 0, 90, 
and 180 degrees using the four phase shifts from Part 3. 


5. Calculate an approximation to the total cross section using these 
four phase shifts. 


Hint: Use the program which initializes the wave function at small r. To 
get better accuracy for the phase shifts, subtract the corresponding 
values of the phase shifts for V = 0. 


Ill. For the potential, V(r) = -exp(-r/.5)/r, and an energy of E = 1: 


rey 


1. Plot the effective potential for 2 = 0, 7, 2, and 3. What is kR for 
this situation? 


2, Plot the radial wave function for 2=O0and X= 7. 
3. Find the phase shifts for 2=0, 1, 2, and 3. 


4. Calculate an approximation to the differential cross section at O, 90, 
and 180 degrees using the four phase shifts from Part 3. 


5. Calculate an approximation to the total cross section using these 
four phase shifts. 


Hint: Use the program which initializes the wave function at small r. To 
get better accuracy for the phase shifts, subtract the corresponding 
values of the phase shifts for V = 0. 


. For the potential, V(r) = 400(1/r?2 - 1/r§), and an energy of E = 1: 


1. Plot the effective potential for 2 = 0, 1, 2, and 3. What is kR for 
this situation? 


2. Plot the radial wave function for 2 = Oand X= 7. 
3. Find the phase shifts for 2= 0,1, 2, and 3. 


4, Calculate an approximation to the differential cross section at O, 90, 
and 180 degrees using the four phase shifts from Part 3. 


5. Calculate an approximation to the total cross section using these 
four phase shifts. 


Hint: Use the program which initializes the wave function at large r. 
For the family of (Mie) potentials, V(r) = 100(1/r™ - 1/r™): 
1, Find the (Q=0, E = 1) and (Q= 1, E = 1) phase shifts for m = 12 and 


n= 1, 2, and 6, Disucss how the behaviors of 5, and 6, reflect the shape 
of the potential. 
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2. Now write the potential as V(r) = 100((rg/r)"* ~ (rg/r)®). Find the 
(2 = 0, E = 1) phase shifts, 5, for rg = 1, 4, 8, and 16, Compare the 
results to 6, from Part 1. Discuss how the behavior of 5 9 reflects the 
characteristic diameter of the potential. 


Hint: Use the program which initializes the wave function at large r. 


. For the family of potentials, V(r) = -exp{-r/rg)/r: 


1. Find the 2 =0, E = 1 phase shift for ro =.1,.5, 1,5, and 10. Discuss 
how the behavior of 6, reflects the screening. 


2. Find the 2 = 1, E = 1 phase shift for ro =.1,.5, 1,5, and 10, Discuss 
how the behavior of 6, reflects the screening. 
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APPENDIX 


INTRODUCTION 


This appendix covers several topics: the terminal plotting performed by the 
subroutines in TTYPLO, the proof that symmetry of a one-dimensional poten- 
tial (for potentials which are finite for finite x) determines that the states can 
be written as wholly even or wholly odd, and a short discussion of the half- 
step method of iterative integration. 


TERMINAL PLOTTING 


The subroutines in TTYPLO (listed below) allow you to plot a set of curves 
on your terminal. The subroutines use the O variables, so you should not use 
these variables in your programs. The output from the subroutines is a 7’’ by 
5" x-y plot (assuming 6 lines/inch and 10 characters/inch) with the x-axis 
across the page and the y-axis up the page. 


You set the range of x values on the plot by defining values for X8 = the left 
side of the plot, X9 = the right side of the plot, Y8 = the bottom of the plot, 
and Y9 = the top of the plot. You then call the subroutines with ““GOSUB 
9000.” 


To plot a point on the graph, you define the values of X@ = the x coordinate 
and Y@ = the y coordinate of the point. You also specify a vatue for ZG, which 
will be the character plotted at the point (X@, Y@) on the graph. You then cail 
“GOSUB 9100.” The allowed values of Z@ are 0 throuth 9. Since these var- 
iables X8, X9, Y8, Y9, X@ Y@, Z@ are used in these special ways for the 
plotting, you should be careful how you use them in your programs. 


To product the final plot on your terminal, you call ‘“‘“GOSUB 9200.” 


A sample of the use of TTYPLO follows. The program plots a sine curve (as 
the character ‘’1’’) and a cosine curve (as the character ‘‘2”’). 


18 READ X8.X9,Y8,Y9 

20 DATA @,6-28318.-1,1 Initialize plot 

38 GOSUB 9960 

48 FOR X@"X8 TO X9 STEP (X9-x8)/70<-—Step across angle 
S6 LET 20*1 

68 LET YO@=SINCX®@) | 
78 GOSUB 918¢ j 
88 LET Z62=2 
96 LET Y@=COS(x®) | Store cosine in plot 
1@@ GOSUB 9100 =} 

ti@ NEXT X@ ; 

128 GosuB 9200<—— Print out plot 


Store sine in plot 


136 STOP 
APPEND-TTYPLO «—— Add p/otting subroutines to program 

Max Y= 1 

XXXKXKXKK KXXKXXKAXXX KX KXKXNRKXKXX KX XK XKKXNKX KKXXKX XK XK NK KKM KX KKM K XK XKXXKKX KXKKX XK 
Y222 1yi101 222 
Y 222 Ht 11 222 
Y 2 lt 11 2 
Y 2 t \ 2 
Y ai ! 2 

Y 12 j 2 

Y i 2 } 2 

Y 1 2 l Z 

oY, ! 2 ] 2 

Y 2 2 

MG 1 2 j 2 

¥. j i 

Y 1 2 1 2 

y 2 2 

Y 1 2 l 2 

Yi 1 } 
% 2 i 2 ! 
% 2 € 

Y 2 1 2 ] 
Y l } 

Y 2 I 2 I 

Y 2 2 

Y 2 l 2 ! 

Y 2 1 2 i 

Y 2 } 2 ] 

Y 2 le l 

Y 2 2h I 

Y 2 a i 

Y 2 2 11 11 

Y 222 222 ti 11 

Y 22222 1bi1i1 


XXKXXXKXKKXXK KX KX KXXXKXK KX XKXAXX KK XXX XAK KX KK X KKK KK XX KKK KK KK XKRKX MK KKKXXAKKAK KK 
MIN Y=~! MIN X® 6 MAX X= 6-28318 
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A LISTing of TTYPLO gives: 


9088S 
92128 
9828 
9836 
9649 
9850 
9968 
9186 
9118 
9128 
9130 
9148 
9188 
9.162 
9172 
9188 
9198 
9200 
9218 
9228 
9238 


9246 
{ 9258 
9268 


9278 
9288 
9298 
9388 
9310 
9328 
9338 
9342 
9352 
9368 
9376 
9388 
9390 
9488 
9416 
9428 
9438 
9440 
9458 
94608 
9a7@ 
9488 
9496 
9586 
9518 
9528 
9538 
9548 
955¢ 
9562 
9570 
9588 
9598 
9 686 
9618 
9628 
9638 
9648 
965@ 
96668 
9678 
9689 


REM INITIALIZE PLOT 
DIM 0€300) 

FOR O!=1 TO 368 

LET 0£ 01392 

NEXT O1 

LET 0€761#76 

RETURN 

REM STORE POINTS 

IF (X8-X8)*"(XB-X9)>8 THEN 9190 . 
IF (Y@-YB)*(YO-Y9)>8 THEN aioe Test for offscale points 
LET 0£763=0f 763+ 1 <-————— Count points plotted 
LET OC OL76)]=1B0G*INT( 3G (YO-¥8) /CY9“Y8) +165) 

LET OC 0€7611 SOC 06 76334 1Q#INT( 704 (XO-X8) /(X9-XB) + 16504722 
IF 0£76)<3@0 THEN 9198~<t———— Jest for too many points 
PRINT "PART OF THE PICTURE NOW" 


GOTO 9200 ~—______Drint out picture 


RETURN 


Initialize the 
storage vector 


Store XO, YO, ZO 
: in one element of 
REM PRINT OUT PLOT 


PRINT "MAX Y="3 ¥9 storage vector 
PRINT ' XXXXXXXXHXXXXXXKMKKXX KAN KX XX XA XK KK AK XXKK KKK XK XAK MAK HKN NK KK X 


FOR 02931 TO 1 STEP -I—«—Step down page 
FOR 03=1 TO 75) 


LET OCO3]=-1 j Initialize line 
NEXT 03 é 
Step through stored points 


LET 0639 
IF 0f033<0 THEN 9369 ——/ est for unused or already used elements 


FOR 03877 TO 0C76) 
LIF INT(0C03)/100@> <> O2 THEN 936¢ «——/s efement in this line? 
LET O1= OC 031~185G*INT< Of 0331/1008) ~«—_ Pamove YG to get XO 


LET OCINTCO1/1@)1201-10"INTCOI/19) 
Remove X@ to get 20 


IF O6>INTCOI/1@) THEN 9352 | 
Cer ctreei | Find largest X@ in line 


LET Of 03) 3-1} qi 
Geen Os Reset element 
RINT "Y's; - 
abi 0381 TO 06 <—olep across line 
GOTO 0£03)42 OF 9488,942859448,9468,9488,9508,9520,9547,9568, 9580 
PRINT " ‘3 
GOTO 9618 
PRINT "'Q@"'3 
GOTO 9618 
PRINT “1"3 
GOTO 9616 
PRINT "2° 
GOTO 9610 
PRINT "3" 
GOTO 96124 
PRINT "4"; 
GOTO 96198 
PRINT "5" 
GOTO 96186 
PRINT ''6°% 
GOTO 9618 
PRINT “7 °'3 
GOTO 9610 
PRINT “8°'3 
GOTO 9618 
PRINT "9 "3 
NEXT 03 
PRINT 
NEXT 062 
PRINT " XXXXXXXXXXKXKXKXXXKXKXX KK KKK KKK RXR KX NK KK KK KKK KK KKK KKK XXX KKK 
"MIN Yu"3YS8s "* MIN X=e"3x8," “os “MAX X8"S X9 


GOSUB 9080 bettas 
“———— Re-initialize plot 


Print correct 
’ character (Z@) 
at each position 


SYMMETRY PROOF FOR ONE-DIMENSIONAL POTENTIALS 


\f the potential is symmetric, then the probability density (the probability per 
unit length), |W|2, will also be symmetric, i.e., |W(-x)|?2 = |W(x)|?. When the 
wave function is real (which usually just means ‘‘when the potential is real’’), 
we have that 


W(-x) = + p(x) (A1) 


as the only possibilities. If Y were a mixture of some even and some odd, then 
lWl2 would not be symmetric. Hence the wave functions can be written as 
wholly even (W{-x) = W(x)) or wholly odd (W(-x) = -W(x)). If there is only 
one wave function for a certain energy, E, then that wave function must be 
either even or odd; it cannot be a mixture. 





A more mathematical proof goes as follows: 


Consider a potential, V, such that V(-x) = V(x). Define an opera- 
tor, P, called the parity operator, by P[f(x)] = f(-x). 


1. The parity operator is Hermitian. 


oo 


LF ¥*60 PIV] dx = f (PL w(x)] * W(x) dy 
where * means complex conjugate. 
2. The eigenvalues of the parity operator are +1 and -1. 
PlW(x)] = a(x) 
P2[h(x)] = aPly(x)} = a2 p(x) 
But P2 [f(x)] = P [f(-x)] = f(x) for any f 
Hence a? = 1 or a = +10r-1 
3. The eigen functions of the parity operator are even and odd. 
From 2: P[W(x)] = W(x) or = Ply(x)] = -W(x) 
for any eigen function 
That means = W(-x) = (x) or = W(-x) = -(x) 


i.e., that the functions are even or odd, respectively. 


4, Finally, then, the parity operator operating on the Schrodinger 
equation with a symmetric potential is the same as the Schro- 
dinger equation for W(-x). 


This in turn implies that the eigen functions of the Schrodinger 
equation can be written as eigen functions of the parity operator. 





THE HALF-STEP ITERATIVE INTEGRATION 


Numerical methods are a subject in mathematics. You may want to study the 
subject someday. For the moment, all we need to show you is that the half- 
step method used heavily in Chapters One and Two would seem to give better 
answers than the simplest approximation (which is called Euler’s method). The 
problem we’re discussing is a general one. A physical law often can be stated 
in terms of the derivative of a function you wish to find; you integrate the 
differential equation to find the answer. 


Consider a general curve for y = f(x); suppose you know the value of y at 
X = Xq and you have a way to calculate the derivative of y with respect to x 
anywhere. The problem is to get as good an approximation to the correct 
value y = f(xg+ax) as reasonably possible. The simplest method (Euler's 
method) uses the derivative at x, and approximates f(x +Ax) as f(x j+Ax) = 
f(xg) + f(x 9) Ax. 
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Euler’s method is shown on the sketch by the dashed line. Clearly, although 
in the limit as x > O the method is correct, for finite Ax the method is not 
very good. 


The geometrical interpretation of the theorem which mathematicians call the 
Mean Value Theorem says that there is some point on the continuously differ- 
entiable curve y = f(x) such that the derivative at that point has the same 
value as the slope of the chord from (x9, f(xg)) to (xg + Ax, f(xq + Ax)). 
if the theorem told us where that point was, everything would be easy. 
So-called higher order methods of iterative integration (such as the fourth- 
order Runge-Kutta method) are ways to approximate better and better this 
“correct’’ value of slope. 


Our half-step method is based on the fact that, in general, evaluating the deri- 
vative near the center of the interval (x9, xg +x) is better than using a deriva- 
tive at the end of the interval. In fact, you can show that, if f(x) is a parabola, 
then a half-step method is exact. This means that, if you expand f(x) around 
Xp in a Taylor series, the half-step method will be correct through terms in 
{Ax)?. That's why half-step methods are called ‘‘second-order’’ Runge-Kutta 
methods, 


wae f(x) 
% step method 








SELECTED ANSWERS 


~ CHAPTER ONE 


Exercise Il. 2(a) E = 2.418 


2(b) none (except t - -%) 
2 2(c) .224 
2(d) 1.41 


3(a) E = 2,507 
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3(b) one (and + - %) 
3(c) 0.389 
3(d) 1.52 


4. E = 089 au. = 2.42 eV = 3.88x10'' ergs = 5.86x10" 
Hertz = 5120A. 


Exercise Ill. 2(a) E = .6845 


fi if 


2(b) none (except + - °%) 
2(c) .812 
2(d) .40 


3(a) E = 2,008 
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3(b) one {and + - °%) 
3(c) 1.294 
3(d) 0.95 


4. E = 1.324 a.u. = 36.03 eV = 5.77x10""" ergs = 8.71x10" 
Hertz = 344A. 


Exercise 1V. 2(a) E =-8.513 





2(b) none (except + - 29) 


2(c) 0.983 
2(d) 0.28 


3(a) E =-5.743 


ie ie 


pL 
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3(b) one {and + - ©) 
3(c) 2.214 
3(d) 0.67 


4. E = 2.76 au. = 75.1 eV = 1.20x10° ergs = 1.82x10'° 
Hertz = 165A. 


Exercise V. 2(a) E = -.407 


2(b) one (and + - ©) 
2(c) 1.407 
2(d) 1.33 


3(a) E =-.0499 





3(b) three (and + - °°) 
3(c) 0.452 
3(d) 6.18 


4. — = .370 au. = 10.1 eV = 1.61x10" ergs = 2.43x107° 
Hertz = 1230 A. 


CHAPTER TWO 


Exercise 1. 2(a) E =-.407 


2(b) E =-.0499 
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2(c} E£ =-.0465 


3. E = .357 a.u. = 9.71 eV = 1.56x10"' ergs = 2.35x10" 
Hertz = 1280A. 


Exercise Hl. 2(a}) E = -.327 


| 
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2(b) E =-.0113 


2(c) E = -.00399 


74 OF QUANTUM MECHANICS 
3 E = .316 a.u. = 8.60 eV = 1.38x10" ergs = 2.08x10"° 
Hertz = 1440A. 


Exercise II!. 2(a) E = -66.27 


J 


2(b) & = -22.98 
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2(c) E = -22.44 


3. E = 43.29 a.u. = 1177 eV = 1.89x10® ergs = 2.85x10"” 
Hertz = 10.5A. 


Exercise IV. 2(a) E = -.0100 


2(b) apparently lost to the continuum 
2(c) apparently lost to the continuum 


3. Since the second state is apparently in the continuum, you 
can set a lower bound on the energy difference. 


E = .0100 au. = .272 eV = 4.36x10''? ergs = 6.58x10" 
Hertz = 45600 A. 
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CHAPTER FOUR 

Exercise 1. 3. Integrating to 20: 
55 = 54.8 - 0.4 = 54.4 degrees 
5, = 42.8 - 1.1 = 41.7 degrees 
6, = 37.7 - 2.4 = 35.3 degrees 
5, = 35.4 - 4.2 = 31.4 degrees 
4. 0 (0) = 4.31 
0 (90) = .0264 
a(180) = .197 
5. O15, = 3.49 

Exercise if. 3. Integrating to 20: 
bg = 124.4 - 0 = 124.4 degrees 
5, = 90.1 - 2.2 = 87.9 degrees 
5. = 74.4 - 6.1 = 68.3 degrees 


5, = 68.2 - 13.0 = 55.2 degrees 


4. o(0) = 91.6 
a (90) = 1.97 

0(180) = 5.98 
5. Oo = 79.9 


Exercise II}. 3. Integrating to 20: 
55 = 29.3 - 0 = 29.3 degrees 
5. = 6.5 - 2.2 = 4.3 degrees 
5, = 7.0 - 6.1 = 0.9 degrees 


63, = 13.2 - 13.0 = 0.2 degrees 


4.0(0) =.317 
0(90) = .104 
0(180) = .058 


5. 0,5, = 1.62 
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